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Abstract. We review recent progresses in the theoretical description of correlation 
C^ ■ and quantum fluctuation phenomena in charge transport through single molecules, 

quantum dots, and quantum wires. A variety of physical phenomena is addressed, 
relating to co-tunneling, pair-tunneling, adiabatic quantum pumping, charge and spin 
fluctuations, and inhomogeneous Luttinger liquids. We review theoretical many-body 
Q ' methods to treat correlation effects, quantum fluctuations, nonequilibrium physics, 

O . and the time evolution into the stationary state of complex nanoelectronic systems. 
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1. Introduction 

In this review we consider different aspects of charge transport through nanoscale 
devices attached to electronic reservoirs, focusing on theoretical approaches dealing 
with interaction and quantum fluctuation effects. Transport experiments have shown 
convincingly that many of these systems - from carbon nanotubes down to single 
nanometer sized molecules - behave as quantum dots. A quantum dot is a confined 
system of electrons, which is so small that the discreteness of energy spectrum with a 
level spacing Ae becomes important; it is therefore often referred to as an artificial atom. 
Due to the smallness of the quantum dot Coulomb interaction between the electrons also 
plays an important role: the charging energy U which can be of same order as the level 
spacing. If the quantum dot is attached to reservoirs by tunneling contacts, electrons 
can leave and enter the dot and the single-particle levels are therefore broadened due 
to the finite lifetime of the electrons. For sufficiently weak tunnel coupling T and small 
temperature T, the level spacing Ae and the charging energy U can be resolved at 
standard cryogenic temperatures. This opens up the possibility of performing transport 
spectroscopy by measuring the differential conductance as function of gate voltage V g 
and bias voltage V. The transport characteristics allow to count the discrete charge 
states of such systems, and even more complex degrees of freedom of the quantum dot 
are revealed. As an example, the spin states can be identified with the help of the 
Zeeman effect, and also vibrational motion of the quantum dot system can be detected. 
The transport behavior, however, depends strongly on the ratios of all the energy 
scales, and on the way transport is induced, e.g., by static or time-dependent applied 
voltages, or by measuring linear or non-linear response in the transport voltage. In 
particular, for molecular scale quantum dots the excitation spectra may be quite complex 
and various splittings appear. This complicates matters further but also leads to a 
variety of transport effects. Besides stationary properties of quantum dots, the time 
evolution out of an initially prepared nonequilibrium state is another important issue 
of recent interest, providing the new energy scale h/t. The time evolution itself can be 
used as a tool to reveal various many-body aspects. It is also of practical importance 
due to the progress in the field of solid-state quantum information processing, initiated 
by the suggestion to use spin states of quantum dots as basic elements of qubits p]. 
In addition to the "bare" energy scales mentioned so far, due to correlation effects 
new effective energy scales can be generated, which describe the physics of quantum 
fluctuations induced by the reservoir-dot coupling. These new energy scales are complex 
functions of the above "bare" energies and we denote them in the following by T c , 
where c stands for "cutoff" . Examples are the Kondo temperature describing the strong 
coupling scale for spin fluctuations (usually denoted by T K ) and renormalised tunneling 
couplings generated by charge fluctuations. These scales will dominate the measured 
transport in the strong coupling regime where T c is the largest energy scale. However, 
also in the weak coupling regime quantum fluctuations set on and are of particular 
interest close to resonances, where they can be analysed by a perturbative treatment 



Charge transport through single molecules, quantum dots, and quantum wires 3 

in the renormalised couplings. To understand the new energy scales physically and 
calculate the renormalised couplings, renormalisation group (RG) ideas are required. 
RG transport theories reveal that interactions are of key importance for the generation 
of new effective energy scales by quantum fluctuations. Without interactions, quantum 
fluctuations are well understood, they just lead to an unrenormalised broadening of the 
energy conservation by the tunneling coupling. Although new effects from correlations 
are most pronounced in the regime of strong correlations r<(7, interesting quantum 
fluctuation effects already start to become visible in the regime of moderate correlations 
r ~ U, where perturbative and mean-field theories are often not yet applicable. Another 
interesting regime is the case of a continuous level spectrum, i.e., where the quantum dot 
is replaced by a one- dimensional quantum wire. In this case, quantum fluctuations in 
the wire itself lead to renormalisation effects and Luttinger liquid physics, with typical 
power-law suppression of the bulk spectral density. Similar to the case of quantum dots, 
the simultaneous presence of correlations and transport through inhomogeneities leads 
to very rich phenomena. A prominent example is the suppression of the spectral density 
at wire boundaries, which leads to vanishing conductance at zero temperature and bias 
voltage if a quantum wire is weakly coupled to two leads. This has to be contrasted to 
the Kondo effect in quantum dots, where conductance becomes maximal when a single 
spin is placed in an "insulating" quantum dot. 

In this review we will illustrate these key issues and concentrate on correlation 
effects in quantum dots and quantum wires. We will describe transport theories 
specifically designed to treat the case of moderate to strong Coulomb interactions 
and quantum fluctuation effects. With these methods it is possible to account for 
the complex interplay of various combinations of quantum mechanical effects (level 
quantization, interference, spin and angular momentum), complex excitations (e.g., due 
to spin-orbit interaction and/or mechanical vibrations), Luttinger liquid physics, non- 
equilibrium (due to static as well as time-dependent applied voltages), and the time 
evolution into the stationary state. An important aspect is that these methods apply 
to general models incorporating experimentally relevant microscopic details, and offer 
prospects for numerical implementation. Despite these common themes, the different 
sections of this review are meant to be almost self-contained and each section closes 
with a short outlook. We consistently use natural units h = k-Q = e = 1. The paper is 
organized as follows. 

In section |2j we will introduce a general approach to transport through quantum 
dots of arbitrary complexity, aiming at the description of three terminal single-molecule 
junctions. The tunneling coupling is treated perturbatively up to second order in 
the coupling T, accounting for charge fluctuation effects, while not capturing spin 
fluctuations which start in higher order. Besides well-known level renormalisation 
and cotunneling effects, electron pair tunneling is shown to give rise to resonances, 
even for the most basic quantum dot model. It is shown how the complexity of, e.g., 
vibrational and magnetic excitations of molecular quantum dots needs to be accounted 
for in transport spectroscopy calculations. 
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Section [3] uses the same perturbation theory but extends it to transport through 
systems which are subject to time-dependent fields. Using the example of adiabatic 
charge pumping through a single- level quantum dot, it explains why time- dependent 
external fields in the adiabatic regime are of particular interest. It is shown how 
time-dependent fields can be used to reveal level renormalisations induced by charge 
fluctuations, an effect purely due to the Coulomb interaction. This offers the interesting 
perspective to use time-dependent fields as fingerprints for the combined influence of 
quantum fluctuations and correlations. 

The regime of strong quantum fluctuations is addressed in section H] in the linear 
response regime. It is described how a recently developed functional RG (fRG) method 
is a useful tool to treat generic quantum dot models. Within this method a systematic 
expansion in a renormalised Coulomb interaction is performed. However, the tunneling 
coupling is treated non-perturbatively, so that the strong coupling regime can be 
addressed for spin as well as for charge fluctuations, at least for moderate Coulomb 
interactions. Two prominent examples of strong quantum fluctuations are reviewed: 
The Kondo effect in quantum dots and the explanation of transmission phase lapses by 
strong broadening and renormalisation effects in multi-level quantum dots. The latter 
involves a significant progress in the insight for an important long-outstanding problem 
posed by experiments. 

Section [5] describes quantum fluctuations in nonlinear response, i.e. at finite bias 
voltage V ^> T, together with the time evolution into the stationary state. The 
basic physics of weak spin and strong charge fluctuations is illustrated on recent 
results for two elementary models: the anisotropic Kondo model at finite magnetic 
field and the interacting resonant level model. A real-time renormalisation group 
(RTRG) method is introduced, where, complementary to the fRG scheme presented 
in section HI an expansion in a renormalised tunneling coupling is performed whereas 
the Coulomb interaction is treated non-perturbatively. This allows the description of 
strong correlation effects r <C U for moderate tunneling. 

Finally, in section |6] we describe transport properties of quantum wires. The status 
of Luttinger liquid physics is reviewed, in particular in connection with the presence of 
inhomogeneities, backscattering, interference effects, and nonequilibrium. As in section 
HI the fRG scheme is shown to be a unique method, capable of describing the whole 
crossover from few- to multi-level dots up to the limit of continuous spectra in quantum 
wires. Its ability to simultaneously incorporate microscopic details is of particular 
importance for the description of experiments. 

2. Transport through single molecule quantum dots 

Molecular quantum dots present the ultimate miniaturization of quantum dot devices, 
which have been realized in semi-conductor heterostructures and wires, metallic 
nanoparticles, and carbon-nanotube molecular wires. Various methods have been 
developed to contact nanometer size single molecules. Transport measurements in 
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three terminal junctions have most clearly demonstrated that single molecules also are 
quantum dots, although of a particularly complex nature, see for a review [2HI]. In view 
of this a central theoretical challenge is to describe the transport through a quantum 
dot of arbitrary complexity, covering the entire, broad class of systems mentioned above. 
We now first outline this general framework, which is relevant also for sections [3] and |5] 
For details see [5] and the references therein. 

Starting from the limit of completely opaque tunnel junction, we can specify the 
exact many-body level spectrum and states \s) of the quantum dot in this case, 

H D =Y,Es\s)(s\. (1) 

s 

Often, only a few of these states are actually accessible in an experiment. Subsequently, 
we incorporate the non-zero transparency of the junction in the coupling 

s,s' aka 

through the amplitudes T^ ka for an electron with spin o to tunnel on / off the dot from 
/ to one of the electronic reservoirs labeled by a = L, R. Generally, in Equation (J2J) all 
states s and s' are coupled which differ by the addition of one electron. The reservoir a 
by itself is described by 

H a = 2_J\ e kv + l^a)n a ka, (3) 

kcr 

where n a ka denotes the electron number operator and k labels the orbital states. Each 
electrode is restricted to contain effectively non-interacting fermion quasi-particles with 
density of states p a . We make the further statistical assumption that the electro- 
chemical potential p a and temperature T are well defined for each reservoir a separately. 
The central quantity is the reduced density operator p of the dot, obtained by partial 
averaging of the full density operator over the reservoirs. From p any non-equilibrium 
expectation value of a local observable can be calculated, and also (see below) the 
transport current. Starting from the above general model, one can derive the kinetic 
equation from which the time-evolution of the dot density operator can be calculated: 

j t p(t) = - iL D pit) -i J dt' £(£, pit'). (4) 

Here the Liouvillian superoperator Ld acts on a density operator by commuting it with 
Hd, i.e., L-£)p(t) = [Hs, p(t)]. This linear transformation generates the "free" time- 
evolution of the dot density matrix, i.e., in the case where it is not coupled to the 
electrodes. The challenge lies in the calculation of the transport kernel £(£,£') which 
incorporates all effects due to the coupling to the electrodes which is adiabatically 
switched on at time to (^o < t). The superoperator £(£,£') is a non-trivial functional 
of three "objects" : (i) the dot Liouvillian L-q (energy differences of initial, final and 
virtual intermediate states of processes), (ii) the dot part of the tunnel operator (tunnel 
vertices) and (iii) the product of the Fermi distribution functions f((eka — fJ>a)/T) and 
density of states p a of all electrodes a. Importantly, compact diagrammatic rules for the 
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Figure 1. (Color online) Sketch of a single-level quantum dot attached to two 
reservoirs with /Ul/r = ±V/2. 



calculation of this kernel in terms of these three objects are known exactly [HIE]. In this 
section and section [3] we review the basic physics which emerges from the calculation of 
the kernel by perturbation theory in the tunnel operator H\ ~ T to leading and next 
to leading order. This is applicable to the "high-temperature" regime, where higher 
order effects are present but still weak. In section [Hand E]the low-temperature physics 
due to strong quantum fluctuations is illustrated. In particular, in section [5] the kernel 
is calculated non-perturbatively in a renormalisation-group framework. The technical 
details of these two approaches can be found in [6] and [5] respectively. 

We now first review all possible basic types of transport resonances which can arise 
due to tunneling processes of the first and second order in T. For this purpose we 
consider a single energy level with Coulomb interaction, the so-called Anderson model, 
where we emphasize the need to consider a finite charging energy U [7] such that double 
occupancy is possible (not only as a virtual intermediate state). This will set the stage 
for the subsequent discussion of examples of transport involving more complex spectra 
in section 12.21 and the time-dependent transport considered in section [3j 

2.1. Transport spectroscopy: sequential, co- and pair-tunneling 

Most types of transport resonances observed in quantum dot experiments can be 
illustrated by the Anderson model sketched in figure [1] where a single orbital level 
is coupled to electrodes a = L, R with tunnel amplitudes T a . The strength of the 
tunnel coupling is characterized by the spectral function, r a = 2irp a \T a \ 2 , whose energy 
dependence is neglected. Denoting by iV the electron number on the dot, the eigenstates 
\s) of the dot (when decoupled from the electrodes) are easily classified. For the empty 
dot (N = 0) there is one state s = with energy E = 0. For the singly occupied dot 
(TV = 1) there are two states s = a =t, \, with energy E a = e a . The N = 1 states 
are split by a magnetic field such that s =\, is the ground state and s =t the excited 
state. Finally, for the doubly occupied dot (N = 2) the energy E& = ^ CT e ff + U of the 
state s = d includes the charging energy. One can show that due to conservation of the 
electron number and the spin-projection along the magnetic field only the occupation 
probabilities enter into the description of the non-equilibrium state. We collect these 
into a vector p = (po,p^,Pi,Pd) T - Its time-evolution is fully described by a kinetic 
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equation with a matrix kernel S (£,£'), 
d 



it p(t) = -i J df'E(t,Op(0 , (5) 

where we have set to = ~~ °° ( c -f- ©)• F° r time-independent parameters considered 
here, we have S(t, £') = £(£ — £') and the occupancies will reach a constant value in the 
stationary long-time limit: p(t) — > p = constant for t — to — ¥ oo. With the subsidiary 
condition of probability conservation Y2 s Ps = 1> this state is determined uniquely by 

= - i S p. (6) 

As expected, here the time integral over the kernel is needed, i.e., the zero- frequency 
Laplace transform lim^io f Q dte lzt H (t, 0) = S. This kernel has been calculated in 
closed form up to order T 2 for an arbitrary dot model Hamiltonian in [6], allowing 
arbitrarily complex transport spectra to be calculated. From a similar equation for the 
measurable tunnel current from electrode a 

I a = -i / dt' Tr D S Q (t, t') p (0 -¥ -i Tr D E Q p, (7) 

J — oo 

the stationary current can be found. This requires an additional zero-frequency kernel 
S a which takes account of both the amount and the direction of the electrons transferred 
from electrode a. Here Trrj denotes the trace over the dot degrees of freedom. 

In figure [2] we show exemplary results for the Anderson model of the calculations 
outlined above. A complete overview of non-linear transport resonances is obtained 
when plotting dl/dV (quantifying changes in the current as new transport processes 
become energetically allowed) as function of the static applied bias V (which drives the 
current) and gate voltage V g (which uniformly shifts all energy levels e CT without changing 
the quantum states) pi In the lower half of the panel we plot the differential conductance 
dl/dV, whereas in the upper half of the panel the occupation p^ of the spin-excited state 
of the quantum dot is shown. This figure illustrates the insight provided by transport 
theory, linking the measured tunnel current to what goes on in the quantum dot. To 
emphasize this correspondence of the two panels, the bias axes in the lower panel is 
mirrored with respect to the upper. 

To the far left and right of the plot the dot is in the N = and N = 2 state 
respectively around zero bias. Moving towards the center from there, dl/dV shows 
a peak along the first encountered yellow skewed lines (marked "SET") where the 
condition fi a = e^ and fi a = e^ + U are met, respectively. Beyond this line a spin 
4- electron can now enter or leave the dot in a sequential or single-electron tunnel (SET) 
process which is of first order in T. In this way the N = 1 ground state s =4 can be 
reached starting from N = and N = 2 respectively. Once beyond the next skew yellow 
line (marked "SET*"), defined by \i a = e^ and fi a = e^ + U, respectively, an electron 
with spin f can tunnel onto the dot. This first order process makes the N = 1 excited 
state s =t accessible, as the upper panel clearly shows. 

| Here one assumes to a good first approximation that the bias and gate electric fields arc spatially 
uniform. Corrections to this picture have been calculated in [8] for molecular junctions. 
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Figure 2. (Color online) Occupation of excited spin-state, p-j- (upper panel) and 
differential conductance (lower panel) for the single-level Anderson model, plotted as 
a function of bias voltage V and gate voltage f3V g , where (3 is the gate coupling factor. 
The spin degeneracy is lifted by an applied magnetic field: e^ — £4. = 50T where T 
is the electron temperature. The dot is symmetrically coupled to the left and right 
electrodes: IY = Tr = 10~ 2 T = 5x 10~ 5 U. The horizontal green line indicates the 
inelastic cotunneling threshold, which equals the Zeeman energy. The skewed green 
lines indicate where the sequential relaxation of the spin-excited state becomes possible 
(COSET). 



Moving further towards the center, there is a large triangular region of the dl/dV 
map where only the N = 1 ground state is accessible: sequential tunneling processes out 
of this state are suppressed to first order since U ^> T, which is referred to as Coulomb 
blockade. However, a second order process, referred to as inelastic cotunneling, allows 
the excited state to become occupied [9j[T0] when the bias exceeds the excitation energy 
for a spin flip, /il — /Ur = V > e-j- — £4.. This defines a gate-voltage independent 
line (horizontal), reflecting the fact that the charge state is left unaltered by the spin- 
fluctuation process which flips the spin. Below this onset voltage only elastic cotunneling 
is possible, giving a smaller conductance and current. In the small triangular region 
marked "ICOT" in the upper panel relaxation of the excited state is possible only by 
inelastic cotunneling. Above this, in the region marked "COSET" , the spin excited state, 
populated by slow inelastic excitation, can relax by a much faster sequential relaxation 
process. It therefore is depleted again and the dl/dV shows a small peak as signature 
of this cotunneling-assisted sequential tunneling (COSET) [TT] . 

The above processes are well-known and have been experimentally observed in 
various types of quantum dots, see, e.g., [9]|12}[T3] and the reviews prfllfll]. However, 
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only recently, it was shown that for this elementary model yet another resonance exists 
in second order in T [7] by taking into account all many-body states and all contributions 
to the transport kernel E. This resonance is related to electron pair tunneling: it shows 
up along a skewed line in the dl/dV map marked "PT", the change in the contrast 
indicating a step in the conductance as function of bias V. The resonance condition for 
this is 2fi a = J2 a £ <? + ^ i indicating that the energy is paid (gained) by two electrons 
coming (going) from reservoir a. This electron pair can be added to (extracted from) 
the dot in one coherent process, thereby changing the charge from N = to N = 2 (or 
vice versa). The electro-chemical potential [i a = \ {Ylo- £cr + ^0 a ^ w hich this happens 
equals the average to the above mentioned SET resonances, i.e., due to quantum charge 
fluctuations the Coulomb energy paid (gained) per electron is reduced by \. 

Finally, we note that the resonance positions of the sequential tunneling peaks are 
slightly shifted with respect to the above mentioned values \i a = e^, e-j- + U due to the 
second order energy level renormalisation which is consistently accounted for. Although 
this may seem unimportant here, in section [3] it will be shown that in time-dependent 
transport this level renormalisation effect may actually generate and dominate the full 
current. 

2.2. Complex excitation spectra: vibrations and spin- orbit splitting 

Molecular quantum dots are clearly more complex than the above considered model, 
for instance by their ability to mechanically vibrate. One important general aspect of 
this complexity is that the only conserved quantity in a tunnel process is the electron 
number |§|. In contrast, the quantum numbers of, e.g., various vibrational modes of 
molecular quantum dot may change when an electron tunnels onto it. This implies 
that the transport theory has to account explicitly for the density matrix off-diagonal 
elements between different excited states with the same charge number N. As is well 
known, this is important already in first order in V when there are degeneracies of the 
energy levels E s in the same charge (N) state on the scale set by the tunneling [T514TT] . 
In the case where there are no such degeneracies only the probabilities matter to order T. 
However, it was shown only recently, that even in this simplest case this is no longer true 
in second (or higher) order in T and Equation ([6]) can no longer be used. Instead, one 
can derive from ([5]) an effective stationary state master equation for the probabilities of 
the same form as ([6]), but with an effective kernel which contains important corrections 
from the off-diagonal elements [6l[T8] . An example where these corrections are found to 
be very important is the Franck-Condon effect. This effect arises in its most elementary 
form in the Anderson-Holstein model, where the dot is described by an Anderson level 
(discussed in section 12.11) plus a single vibrational mode with frequency to: 

H D = ^ £ff n (T + £/ nt% + Av / 2Q^ CT + !(P 2 + Q 2 ). (8) 

cr a 

§ Conserved here refers to the system of reservoirs and dot, including their interaction 
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Figure 3. (Color online) Differential conductance map vs. applied bias V and gate 
voltage V g {fi — gate coupling factor) for a molecular level coupled to a vibrational 
mode with frequency u> and strong local interaction U 3> T, V, A [B] . With increasing 
electron vibration coupling A, the inelastic cotunncling steps (red arrows in the A = 2 
figure indicate 1 and 2 phonon absorption) become relatively important with respect 
to the suppressed SET processes. The white lines / regions correspond to negative 
dl/dV, which cannot be plotted due to the logarithmic color scale. Strikingly, the 
current thus displays a peak in the Coulomb blockade regime, which is very uncommon. 
In all plots we have chosen u> = 40T = 10 4 r and conduction bandwidth D = 250a;. 
The conductance is scaled to Tm, the maximal sequential tunneling rate, i.e., T times 
the maximal Franck-Condon overlap factor (which is less than 1). 



Here Q = (tf + b) /a/2 is the vibrational coordinate normalized to the zero-point motion 
and P = i(tV — 6)/a/2 the conjugate momentum. Due to the linear coupling term oc A 
the vibrational mode is displaced when charging the molecule. For A > 1 this causes 
a suppression of the Franck-Condon overlap of ground-state vibrational wave functions 
for two subsequent molecular charges N. As a result the tunnel rates calculated to 
lowest order in T are suppressed and the sequential tunneling current is blocked, even 
for gate voltages close to the resonance [T9H2T| . an effect also called "Franck-Condon 
blockade" [22]. This implies that second order processes become of crucial importance 
even at resonance [23J. The general approach developed in [6] consistently accounts 
for all second order effects, i.e., not just cotunneling rates as in [23]. In figure [3] 
we show the differential conductance calculated to second order where the vibrational 
mode is shifted by several multiples of its zero-point amplitude. Whereas for A = 1 the 
standard Coulomb blockade picture is still identifiable, for A = 4 it is drastically altered: 
inelastic cotunneling (threshold at /xl — A*r = w, 2w, . . .) and COSET [23] resonances 
proliferate in the transport spectrum. Experimentally, vibration assisted transport in 
three terminal molecular devices has been reported [251427] . Significant Franck-Condon 
blockade has been observed in suspended carbon nanotubes and compared with theory 
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in [22]. Other theoretical studies have studied the effect on current noise [221120], the 
Kondo effect [301131], the influence of mechanical dissipation (finite Q-factor) [32],[33] or 
the lack thereof [M], vibrational anharmonicity [35j[36] and interference effects [351137] . 

All above works rely on the Born-Oppenheimer separation of nuclear and electronic 
motion on the molecule. Transport effects signaling the breakdown of this separation, 
going under the generic name of (pseudo-) Jahn- Teller dynamics, have also been 
calculated [38H1U] . A very recent STM experiment jH] has confirmed these transport 
signatures of the breakdown of the Born-Oppenheimer approximation. In this work SET 
resonances were measured as function of a molecular parameter, as suggested in [3D], in 
this case by measuring many individual molecular wires of various lengths. 

Another class of examples of quantum dots with complex excitations are magnetic 
molecules. Here the spin-orbit interaction results in violation of spin-selection rules for 
electron tunneling and the same remarks as above apply. In particular, single molecule 
magnets, reviewed in [12], offer a basic example and have been of interest in view of 
quantum computing P2HHJ. In such molecules the ground state has a large spin value 
S > 1/2 due to strong exchange coupling of the spins of a few magnetic atoms. This 
multiplet is however split by strong spin-orbit perturbations. In many cases an accurate 
model for a given charge state N is 

<> = - D N Sl + E N {Si - Si) , (9) 

where the uni-axial (Dn) and transverse (-EW) anisotropy parameters describe the 
splitting and mixing of spin states, respectively. Recently, three-terminal measurements 
on magnetic molecules have been reported, including comparison with theory [45-48J, 
see for a review [10]. Calculations of the sequential [5DTI55] and cotunneling [5&H6D] 
regimes for models including the anisotropy have identified distinctive transport 
signatures and possibly useful effects. An important goal is to gain control over 
the spin of a single molecule by electric means [611 [62]. Recently, the magnetic 
anisotropy splittings were measured in a three terminal junction, similar to two terminal 
STM measurements [631 El]- Importantly, the magnetic field evolution of inelastic 
cotunneling [18] was measured in different charge states. It was found that the magnetic 
Dn parameter depends strongly on the charge state N, which can be electrically 
controlled with the gate voltage. 

Beyond the issue of control over single-molecule magnetism, is the understanding 
of coupled molecular magnets or magnetic centers. Apart from spin-blockade effects, 
expected for large spin molecular ground states [65], the effect of spin-orbit interaction 
on the exchange coupling between two magnetic centers inside a magnetic double 
quantum dot is of interest. For instance, the spin-orbit induced Dzyaloshinskii-Moriya or 
antisymmetric exchange interaction was shown to give rise to a characteristic dependence 
of the transport on an externally applied magnetic field and the polarizations of magnetic 
electrodes [66] . 

Higher order tunneling effects for magnetically anisotropic molecular quantum 
dots described by fl9]) have attracted attention [6TH69] after the prediction of a novel 
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anisotropic pseudo spin Kondo effect [70]. The spin-orbit induced anisotropy in Q not 
only raises a barrier which opposes spin-fluctuations (due to D N ) (section 0J, but it 
also generates intrinsic spin-tunneling (due to En). This interplay leads to interesting 
anisotropic effective exchange interaction with electrons transported through a molecular 
magnet, as will be discussed in section [5] 

Besides their own interesting transport signatures, spin and vibrations can also 
interact or influence each other in transport. For instance, a mechanism of vibration- 
induced spin-blockade of transport was predicted for a mixed-valence dimer molecular 
transistor ("double dot") [7T|I72]. 

Finally, we mention that level spectrum complexity may also deceive one when 
analyzing experimental data. A particularly clear example is offered by recent transport 
data on carbon-nanotube "peapods", i.e. fullerene molecules in a host nanotube |73j . 
Here weakly gate-dependent dl/dV peaks were observed which are strongly reminiscent 
of inelastic cotunneling (second order in T, see section |2~T|) . However, detailed transport 
calculations have conclusively shown that the complex measured spectrum can be 
assigned to sequential tunneling processes only, i.e., of first order in T. This deceptive 
case occurs because the two coupled quantum dots, the host nanotube and the fullerene 
"impurities" have a very different gate voltage dependence. Fortunately, the significant 
hybridization of the two dots gives a characteristic interference effect in the first order 
current, distinguishing it clearly form second order effects. In view of applications it 
is interesting that the analysis in [73] indicates that the charge state of the fullerene 
impurities is electrically tunable, independently of that of the host carbon nanotube. 

2.3. Outlook 

In summary, we have illustrated that general methods for spectroscopy calculations 
for molecular quantum dot models are essential for the experimental characterization 
and ultimately the design of nanoscale quantum dot devices. The complex spectra 
due to various strong many-body interactions play a crucial but complicating role, 
at the same time however, enabling interesting new possibilities. An exciting future 
direction important for further progress is to efficiently account for renormalisation 
effects in complex models, in particular as outlined in section [5] Also, the approach 
outlined here can be extended to other transport quantities, such as spin and heat [71] . 
The molecule and junction models used here can be parametrized based on atomistic 
electronic structure calculations, as for instance, in [I7J[75j[76] . Pursuing this further is 
another important step in understanding molecular quantum dots. 

3. Adiabatic transport through a quantum dot due to time-dependent fields 

Charge transport through quantum dots is not necessarily due to an applied static bias 
voltage: in the absence of a bias, a current can be obtained by the time-dependent 
modulation of fields, applied externally to a mesoscopic device. This is of interest for 



Charge transport through single molecules, quantum dots, and quantum wires 13 

possible future applications in nanoelectronics, where the repeated and fast operation 
of a device is of interest. Indeed the mesoscopic analogue of a capacitor has been 
realized [77J,[7E] serving as a fast coherent electron source and charge pumping [791482] is 
a candidate for the realization of a quantum current standard, to name some examples. 
From a fundamental point of view the nonequilibrium created by time-dependent electric 
fields provides enhanced insight into the quantum properties of a nanoscale device, which 
are not accessible with equilibrium methods. If the variation of the fields is slow with 
respect to the internal time-scales of the system itself, the modulation is called adiabatic. 
The system then remains in a quasi-equilibrium state, meaningthat the system is not 
brought to an excited state by the time-dependent modulation l]j]. 

A prominent example for transport in the absence of a bias voltage is adiabatic 
pumping: the cyclic variation of at least two of the system parameters, see figure HJ 
leads to a dc charge transfer through the quantum device. In this case the pumped 
charge does surprisingly not depend on the detailed time evolution of the pumping 
cycle; this can be shown by formulating it in terms of geometric phases [83H85] . 
Theoretically, adiabatic pumping has been extensively studied in the non-interacting 
regime [861 - 192"] . In this case a scattering matrix approach is convenient j86j|93] and has 
been broadly used. While indeed Coulomb interactions are of minor importance in some 
setups [TTJEHlini], they play an important role in many quantum-dot systems and lead 
to a variety of interesting results. Several approaches have been used to address the 
problem of adiabatic pumping through interacting quantum dots in different regimes, 
including weak Coulomb interaction j95j|96], pumping in Luttinger liquids [97], and 
a quite general approach using nonequilibrium Green's functions [981 - 1100] . The latter 
works are particularly appropriate for the study of weak interaction or for the calculation 
of spin and charge pumping in the Kondo regime [99] and the role of elastic and inelastic 
scattering processes has been studied in detail [100] . The treatment of the tunnel 
coupling is non-perturbative in these works and they are in this sense complementary 
to the results presented in the following. 

It turns out that the interplay of Coulomb interaction and coherent tunneling 
between a quantum dot and leads, leading to an effective renormalisation of the energy 
level, can be at the origin of the pumping mechanism. The result is that this energy 
renormalisation, which is a pure Coulomb interaction effect vanishing when Coulomb 
interactions are negligible, is accessible via the pumped charge. 

Here we present results in the regime in which the coupling between the dot and 
the leads is moderate; this means that the energy scale T approaches temperature T 
but is still smaller than T so that charge fluctuations start to become important. In 
contrast, the energy scale T K of the Kondo temperature setting the scale for the onset 
of spin fluctuations, is taken to be much smaller than T and T. To approach the effect 
described above we use a perturbative expansion in the tunnel coupling between the 

j| The slow variation of the externally applied fields becomes experimentally particularly important 
for delicate systems as molecules, see section [2 where problems related to heating, surface excitations 
(plasmons) and laser-induced chemical reactions (photo-bleaching) should be avoided. 
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Figure 4. (Color online) (a) Sketch of a single-level quantum dot with time-dependent 
parameters attached to two reservoirs, (b) Two of the system parameters are time 
dependent as indicated in (a) and perform a closed cycle in parameter space. 



dot and the leads up to second order in T, which is appropriate in this regime. The 
aim is the understanding of the influence of finite and arbitrary Coulomb interaction on 
the pumping characteristics and the identification of the physical origin of the various 
contributions to the pumped charge. In order to achieve this, a diagrammatic real- 
time technique [101H104] . which has been developed to describe non-equilibrium DC 
transport through an interacting quantum dot, see section (J2J), was extended for the 
treatment of adiabatically time-dependent fields |105] . 



3.1. Adiabatic expansion 

Even though the approach is not restricted to specific systems, here the focus is put on a 
single-level quantum dot as shown in figure HJ The quantum dot has a time-dependent 
level e(t) = e + Se(t) and is coupled to leads a = L, R with time-dependent tunnel 
amplitudes T a (t). The strength of the tunnel coupling is characterized by the time- 
dependent quantity, T a (t) = 27rp a T*(t)T a (t) = f a + ST a (t), with the lead's constant 
density of states, p a . As described in section [2j the eigenstates of the decoupled 
system are given by s = {0, t, \., d}. Tracing out the lead degrees of freedom the 
system is described by its reduced density matrix, where here only the occupation 
probabilities p = {po,p^,Pi,Pd) T axe of interest. Their time-evolution in the long-time 
limit, to —* — oo, is fully described by the generalized Master equation, Equation (JSJ. 
The goal being the description of slowly varying system parameters, it is useful to 
perform an adiabatic expansion |105j : the frequency of the variation therefore has to 
be much smaller than the dwell time of electrons in the system, 0< T. The effect of 
the system parameters varying slowly in time is that the density matrix of the system 
at time t is not only described by the instantaneous values of the parameters, but it 
lags a bit behind the parameter change. Therefore an expansion around this time t 
suggests itself. To obtain such an expansion as a first step a Taylor expansion of p(t') 
is performed around t up to linear order in the integrand on the right hand side of 
Equation (jHj). This is related to the fact that memory effects of the kernel have to be 
taken into account. Furthermore, an adiabatic expansion of the kernel S (£,£') itself is 
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Figure 5. (Color online) (a) Pumped charge due to a parameter variation of the dot 
level and of one of the barriers as a function of the average dot level. The dominant 
contribution to the pumped charge is due to first order tunneling processes, (b) 
Pumped charge due to the pure variation of the barrier strengths as a function of 
the average dot level. The dominant contribution to the pumped charge is due to 
second order tunneling processes only. In both plots the temperature equals T = 2T. 



performed. The zeroth-order term, Sj [t — t'), is indicated with the superscript (i) for 
instantaneous. The subscript t to emphasize that the system parameters X{j) — > X(t) 
are frozen at time t, i.e., the functional dependence on X(r) is replaced by X{r) — > X(t). 
The first-order term is obtained by linearizing the time dependence of all parameters 
X(t) with respect to the final time t, i.e., X(r) — > X(t) + (r — t)-^X(r)\ T=t , and 
retaining only linear terms in time derivatives. This linear correction to the kernel is 
indicated by the superscript (a) for adiabatic. Finally, it is necessary to perform an 
adiabatic expansion for the occupation probabilities in the dot, 
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The instantaneous probabilities p\ are the solution of the time-independent problem 
with all parameter values fixed at time t. They are obtained by solving the stationary 
Master equation, Equation (jSj), with all parameters replaced by time-dependent 
parameters evaluated at time t. Once the instantaneous probabilities pf are known, 
the adiabatic corrections pf are found from the Master equation in first order in Vt 
obtained from the expansion described above. Similarly an equation for the current can 
be developed, see equation (J7|), and adiabatically expanded, where a current kernel S L 
takes account for the amount and the direction of electron transfer. 

On top of the adiabatic expansion a perturbative expansion in the tunnel coupling 
strength V is performed up to second order, allowing the actual evaluation of the 
kernel using diagrammatic rules as elaborated in Refs. |101H105] , see section [2j This 
perturbative expansion is valid in the regime of moderate couplings with respect to 
temperature, i.e. Y <T. 
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3.2. Pumping current 

We are interested in the time-resolved pumping current as well as the charge pumped 
through the dot per pumping cycle. The pumped charge, Q = L Ii,(t)dt, in the bilinear 
response regime, i.e. for modulations with small amplitudes, is proportional to the area 
spanned in parameter space, see figure HI As a first step the time-resolved pumping 
current is calculated in lowest order in the tunnel coupling. It is straightforward to 
relate the pumped current to the dynamics of the average instantaneous charge of the 
dot (n)t and one finds 

li 0) (t) = -^##(n)f' 0) , (12) 

L v ; r(t) dv n ' v ; 

where {n)t denotes the instantaneous occupation in lowest order in T. This suggests 
the following interpretation. As the dot occupation in lowest order, (ra)f' , is changed in 
time by varying the pumping parameters the charge moves in and out of the quantum dot 
generating a current from/into the leads. Note that one of the pumping parameters must 
be the level position since (n)i h is given by Boltzmann distributions and is therefore 
independent of the tunnel-coupling strengths. This means that ^(n)t = if the level 
position is constant. The charge moving in and out of the quantum dot is divided to 
the left and to the right, depending on the time-dependent relative tunnel couplings 
T a (t)/T(t). 

The next step is to calculate the next-to-leading order correction to the lowest 
order result shown in Equation ( JT2l) . Also the next-to-leading order can be written in 
a compact way 

r (l) U\ - I ( I \(*,broad,L)\ ^(t) d . , (i,ren) 



^W = -1SI<">' -^J + TfiW' V (13) 

The first term of Equation (fl3j) contains the contribution due to the correction of 
the average dot occupation induced by lifetime-broadening. It contains a total time 
derivative, and as parameters are periodically changing in time, it will not lead to a net 
pumped charge after the full pumping cycle. The second term has the same structure 
as the zeroth-order contribution, Equation (IT2"|) . It can be understood as the correction 
term induced by the combined influence of charge fluctuations and correlation effects 
giving rise to a renormalisation of the level position, e{t) — > e(t) + a (e(t),T(t),U). 
The level renormalisation, a (e,T,U), is positive when the level is in the vicinity of 
the Fermi energy of the leads and negative for e + U being close to the Fermi energy. 
The sign of the renormalisation thus indicates if transitions take place between empty 
and singly occupied dot or between single and doubly occupied dot. The energy level 
renormalisation may be time dependent via time-dependent tunnel couplings or a time- 
dependent level. The result is that a finite charge can be pumped by means of level 
renormalisation. This is a pure Coulomb interaction effect, it vanishes for U — Y 0. We 
note that correction terms from the renormalisation of the tunneling couplings do not 
contribute in first order in V. For the DC current, different contributions from higher 
order processes are present at the same time, which makes it challenging to identify them 
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separately in an experiment. For the adiabatically pumped charge, where correction 
terms associated with a renormalisation of the tunnel couplings and level-broadening 
effects vanish, the situation is distinctively different. Studying adiabatic pumping is, 
therefore, a convenient tool to access the energy-level renormalisation. This is most 
dramatic in the case when the zeroth-order pumped current is zero, i.e., when pumping 
is done by changing both tunnel couplings. In this case, the dominant contribution 
to the pumped charge is due to time- dependent level renormalisation. These results 
suggest that adiabatic pumping can be used to directly access the level renormalisation 
in quantum dots. 

The results for the pumped charge per pumping cycle are shown in figure El where rj 
is the enclosed surface in parameter space. Modulating the level position and one of the 
tunnel couplings the pumped charge has a maximum contribution at the resonances, 
see figure [5(a). Importantly the contribution from sequential tunneling is dominant 
here. Figure [5(b) shows the pumped charge obtained by modulating the two tunnel 
couplings. As described above, the pumped charge is due to level renormalisation and 
therefore vanishes for vanishing Coulomb interaction. The sign change between the two 
contributions at the resonances reflects the opposite sign of the level renormalisation for 
the two resonances. 

The method described here is generally formulated and it is therefore extendable 
to a variety of systems in which Coulomb interaction is important and to which time- 
dependent fields are applied. Charge and spin pumping through an interacting quantum 
dot has been studied in the presence of ferromagnetic leads |106j . Pumping through 
one or two metallic interacting islands with a continuous density of states has been 
examined in |107j . In |108] the validity of the frequency regime has been extended to 
faster modulations and in |109j the influence of quantum interference was studied on 
the pumped charge through one or two quantum dots embedded in an Aharonov-Bohm 
geometry. Finally this approach has been used to study the capacitive and the relaxation 
properties of a driven quantum dot figuring as a mesoscopic capacitor in |110j . 

3.3. Outlook 

In this section we have presented results on adiabatic pumping, where interesting effects 
due to quantum charge fluctuations and finite Coulomb interaction are revealed. The 
underlying method shown here uses an adiabatic expansion of the Master equation based 
on a real-time diagrammatic technique; this method is applicable whenever a mesoscopic 
system is exposed to a certain number of adiabatically time-dependent fields. The 
intriguing results found for systems with a relatively simple spectrum together with the 
general formulation of the adiabatic expansion of the Master equation, motivate further 
investigations. A challenging question is, e.g., to study the impact of time-dependent 
fields on molecular devices with a complex energy spectrum, as introduced in section [2j 
Another interesting development concerns the combination of the formalism of adiabatic 
quantum pumping with renormalisation group methods, as described in sections H] and 
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to describe the influence of time-dependent external fields in the regime of strong 
correlations and strong quantum fluctuations. 

4. Quantum fluctuations in linear response 

In this section we will discuss the physics of strong quantum fluctuations in combination 
with correlation effects in quantum dots. We will concentrate on the linear response and 
static regime, the dependence on finite bias voltage together with the time evolution 
will be discussed in section [5j One of the most prominent examples of the drastic 
effects of spin fluctuations in quantum dots is the experimental observation of the 
Kondo effect [111] . In bulk solids a small amount of magnetic impurities leads to 
an increased magnetic scattering of the electrons at low temperature, which results 
in an increased resistance. In quantum dots, the mesoscopic realization of a single 
spin coupled to two leads displays instead a zero-bias peak in the conductance |112j . 
The particular advantage of quantum dots is that they allow for a very flexible tuning 
of the parameters and can easily be extended to study the effects of more complex 
impurities, where orbital Kondo, quantum critical and interference effects may arise. 
Moreover, additional environmental degrees of freedom in presence of ferromagnetic 
or superconducting reservoirs coupled to the quantum dot as well as finite-size effects 
affect the transport properties. For spin fluctuations, the characteristic energy scale T c 
is given by the Kondo temperature Tk, which, for the special case of a single- level dot 
with Coulomb energy U and tunneling coupling T, is given by Tk ~ yUYe~^c . This 
scale is exponentially small for large Coulomb energy and, therefore, the Kondo effect in 
quantum dots is only visible for strongly coupled leads. In contrast, the importance of 
charge fluctuations is controlled by V and signatures are already significant for T ~ T. 
Already in the last section, we showed how renormalised level positions can be identified 
in the adiabatically pumped charge. Experimentally, effects from charge fluctuations 
have been detected in transmission phase lapses through multi-level quantum dots in 
Aharonov-Bohm geometries |113[I114] . New light was shed on this long-outstanding 
puzzle by the insight that the combined influence of broadening and renormalisation 
effects induced by charge fluctuations is responsible for this effect jl 15(,lll6j . 

The analysis of the signatures of correlation and quantum fluctuation effects 
requires appropriate methods for their theoretical description. The importance of the 
development of analytical as well as numerical techniques for their treatment constitutes 
therefore an important issue. In recent years functional renormalisat ion-group (fRG) 
methods [117] have been established as a new computational tool in the theory of 
interacting Fermi systems. These methods are particularly powerful in low dimensions. 
The low-energy behavior usually described by an effective field theory can be computed 
for a concrete microscopic model by solving a differential flow with the energy scale 
as the flow parameter. Thereby also the nonuniversal behavior at intermediate scales 
is obtained. For applications to quantum dots and wires the fRG scheme turned out 
to be an efficient approach for the description of the single-particle spectral properties 
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and the linear transport through generic models in presence of moderate correlations. 
Arbitrary tunneling couplings can be considered, allowing to access the regime of strong 
coupling between the reservoirs and the quantum system. Due to the flexibility of 
the microscopic modeling the whole parameter range can be easily explored, and more 
complex geometries of multi-level quantum dots, as realized in recent experiments, can 
be treated. In the following sections we will discuss applications of the fRG to the 
description of strong spin and charge fluctuations in quantum dots, exemplified by the 
discussion of the Kondo effect and transmission phase lapses in multi-level quantum 
dots. 

4-1. The Kondo effect in quantum dots 

The appearance of an upturn followed by a low-temperature saturation of the resistivity 
in metals containing diluted magnetic impurities was explained by Kondo in 1960 as an 
enhanced scattering due to the screening of the local impurity spins by the conduction 
electrons |118j . The Fermi liquid nature of the ground state |119] and the various 
theoretical treatments, including Wilson's numerical renormalisation group (NRG) [120j . 
make the Kondo problem one of the best understood many-body phenomenon in 
condensed matter physics, as well as a paradigm for correlated electron physics |121j . 
In the late 90's, the prediction by Ng and Lee |122j . as well as Glazman and Raikh |112] 
of the Kondo effect in quantum dots led to a revival of the subject, and has been 
observed in 2D GaAs/AlGaAs heterostructure quantum dots |llllll23lll24| . silicon 
MOSFET's, as well as in carbon nanotubes |125] and contacted single molecules 
[126(1127] . Confinement and electrostatic gating enabled to investigate the Kondo 
effect in an artificial Anderson impurity through electrical transport measurements, 
characterized by a unitary conductance 2 y at zero temperature |128] . Besides providing 
a paradigm for a variety of physical effects involving strong electronic correlations, the 
Kondo effect in nanostructures allows for the realization of spintronic device constituents 
as spin filters or quantum gates. The understanding and control of the transport of spin- 
polarized currents are of fundamental importance for the advancement of semiconductor 
spintronic device technology |129j . 

The coupling of a quantum dot with spin-degenerate levels and local Coulomb 
interaction to metallic leads gives rise to Kondo physics |121j . A small quantum dot with 
large level spacing is described by the single-impurity Anderson model (SIAM) in the 
regime in which only a single spin-degenerate level is relevant. At low temperatures and 
for sufficiently high barriers the local Coulomb repulsion U leads to a broad resonance 
plateau in the linear conductance G of such a setup as a function of a gate voltage V^ 
which linearly shifts the level positions |112[I122[I1301 - I132"] . It replaces the Lorentzian of 
width T found for noninteracting electrons. On resonance the dot is half-filled implying 
a local spin-i degree of freedom responsible for the Kondo effect |121j . In the limit of 
large U ^> T the charge degrees of freedom can be integrated out and the spin physics is 
described by the Kondo model, which will be discussed in section [5j For the SIAM the 
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zero temperature conductance is proportional to the one-particle spectral weight of the 
dot at the chemical potential [133] . The appearance of the plateau of width U in the 
conductance is due to the pinning of the Kondo resonance in the spectral function at the 
chemical potential for — 2 < V g < ^ (here Vg = e + ^ = corresponds to the half-filled 
dot case) [121U132J . Kondo physics in transport through quantum dots was confirmed 
experimentally jlllj.1134] . and theoretically using the Bethe ansatz |130[I132] and the 
NRG technique [12011135] . However, both methods can hardly be used to study more 
complex setups. In particular, the extension of the NRG to more complex geometries 
beyond single- level quantum-dot systems |131[I132[I136] is restricted by the increasing 
computational effort with the number of interacting degrees of freedom. Alternative 
methods which allow for a systematic investigation are therefore required. Here the fRG 
approach is proposed to study low-temperature transport properties through mesoscopic 
systems with local Coulomb correlations. 

A particular challenge in the description of quantum dots is their distinct behavior 
on different energy scales, and the appearance of collective phenomena at new energy 
scales not manifest in the underlying microscopic model. An example of this is the 
Kondo effect where the interplay of the localized electron spin on the dot and the spins 
of the lead electrons leads to an exponentially small (in U/T) scale Tk- This diversity 
of scales cannot be captured by straightforward perturbation theory. One tool to cope 
with such systems is the renormalisation group: by treating different energy scales 
successively, one can often find an efficient description for each one. The fermionic 
fRG |137] is formulated in terms of an exact hierarchy of coupled flow equations for 
the vertex functions (effective multi-particle interactions) as the energy scale is lowered. 
The flow starts directly from the microscopic model, thus including nonuniversal effects 
on higher energy scales from the outset, in contrast to effective field theories capturing 
only the asymptotic behavior. As the cutoff scale is lowered, fluctuations at lower energy 
scales are successively included in the determination of the effective correlation functions. 
This allows to control infrared singularities and competing instabilities in an unbiased 
way. Truncations of the flow equation hierarchy and suitable parametrizations of the 
frequency and momentum dependence of the vertex functions lead to new approximation 
schemes, which are devised for moderate renormalised interactions. A comparison with 
exact results shows that the fRG is remarkably accurate even for sizeable interactions. 

To exemplify the effect of correlations, we first focus on a quantum dot with spin- 
degenerate levels. For simplicity only a single level with a local Coulomb repulsion U 
described by the SIAM |121] is considered, see figure [TJ The hybridization of the dot 
with the leads broadens the levels on the dot by r a = 2-irT^p a , where p a is the local 
density of states at the end of lead a = L,R assumed to be independent of frequency 
here. The energy level is determined by the gate voltage V g . Integrating out the lead 
degrees of freedom, the bare Green function of the dot |121j is 

SoM = .. \ r : : , (14) 

ibj - Vg + 2|sgn(w) 

where r = i\ + Tr. By solving the interacting many-body problem a self-energy 
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Figure 6. (Color online) Upper panel: conductance as a function of gate voltage for 
different values of U/T. Lower panel: average number of electrons on the dot. 



contribution T,(iu) is obtained dressing the bare propagator through the Dyson equation. 
The fRG is used to provide a self-energy describing the physical properties of the T = 
linear conductance obtained by G(V g ) = ^-7iTp(0) |133] in terms of the dot spectral 
function p{u) = — -ImQ(cu + i0 + ). In order to implement the fRG, an infrared cutoff 
in the bare propagator is introduced Qq(zu) = G (iu))Q(\tu\ — A). As the cutoff scale A 
is gradually lowered, more and more low-energy degrees of freedom are included, until 
finally the original model is recovered for A — > 0. Changing the cutoff scale leads to an 
infinite hierarchy of flow equations for the vertex functions. In the static approximation 
the flow equation for the effective level position V A = V g + S A reads 

A _ U A V A /n 



d A V 



(15) 



(A + 2 + (K A ) 2 

with the initial condition V A= °° = V g |138j . In first approximation the two-particle 
vertex U A = U A=OD = U equals the bare Coulomb interaction. At the end of the flow, 

rA=0 



the renormalised potential V = V determines the conductance G(V g ) = 2 ^- 4V2+T 2 ■ 
Results for different values of U/T are shown in figure El together with the occupation 
of the dot. 

For r< (7 the resonance exhibits a plateau |132j . In this region the occupation 
is close to 1 while it sharply rises/drops to 2/0 to the left/right of the plateau. For 
asymmetric barriers the resonance height is reduced to 2 ^- (r _P r R ^ 1121(1132] . Focusing 

1U 

on strong couplings U 3> T, the solution of the above flow equation V ~ V^e~^r 
describes the exponential pinning of the spectral weight at the chemical potential for 
small \V g \ and the sharp crossover for a V of order U. While already the first order 
in the flow-equation hierarchy captures the correct physical behavior, the inclusion of 
the renormalisation of the two-particle vertex improves the quantitative accuracy of 
the results. For details on the parametrization and extensions including dynamical 
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properties we refer to Refs. [138H142] . In presence of a magnetic field the Kondo 
resonance splits into two peaks with a dip in G(V S ) at V g = 0, providing a definition of 
the Kondo scale as the magnetic field required to suppress the total conductance to one 
half of the unitary limit. For the single dot at T = the conductance and transmission 
phase are related by a generalized Friedel sum rule to the dot occupancy [121] by 
G = 2 2- sin 2 (| (n)) and a — f (n). For gate voltages within the conductance plateau 
the dot filling is 1 and the phase is |. The description of transmission phases for more 
complex setups will be discussed in the next subsection. 

4-2. Mesoscopic to universal crossover of the transmission phase of multi-level 
quantum dots 

The following application to a multi-level quantum dot illustrates the strength of the fRG 
approach, the flexibility and simple implementation, in the description of the intriguing 
phase-lapse behavior observed in experiments by the group of Heiblum at the Weizman 
Institute |113]I114] . The transmission amplitude and phase T = |T|e iQ of electrons 
passing through a quantum dot embedded in an Aharonov-Bohm geometry showed 
a series of peaks in \T\ as a function of a plunger gate voltage V g shifting the dot's 
single-particle energy levels. Across these Coulomb blockade peaks a(V g ) continuously 
increased by it, as expected for Breit-Wigner resonances. In the valleys the transmission 
phase revealed "universal" jumps by tt for large dot occupation numbers. In contrast, 
for small dot fillings the appearance of a phase lapse depends on the dot parameters in 
the "mesoscopic" regime. From the theoretical side, the observed behavior is captured 
by an fRG computation of the transport through a multi-level quantum dot with local 
Coulomb interactions [115(1116] . An essential aspect in the description of the generic 
gate- voltage dependence of the transmission lies in the feasibility of a systematic analysis 
of the whole parameter space. The interaction is taken into account approximately, but 
a comparison to numerically exact NRG data for special parameters proves the fRG 
results to be reliable as long as the interaction parameter and the number of almost 
degenerate levels do not become too large simultaneously. The results are shown in 
figure 0, for spin-degenerate levels we refer to Ref . [116] . 

Assuming level spacings of the order of the level width or smaller for large dot 
fillings (similar to a hydrogen atom with decreasing level spacing for increasing quantum 
number) and well separated levels for dots occupied by only a few electrons, the 
obtained results are consistent with the experimental findings in both regimes. In 
the "mesoscopic" regime the Coulomb blockade yields an increase in the separation 
of the transmission peaks. The behavior in the valleys depends on the details of the 
level-lead couplings, and can be continuous or discontinuous with a phase lapse of n 
(see upper panels in figure [7]). For several almost degenerate levels in the "universal" 
regime the hybridization leads to a single broad and several narrow effective levels (Dicke 
effect [143] ). In presence of the Coulomb repulsion, the gate- voltage dependence of the 
broad level is significantly reduced, in contrast to the narrow levels crossing the broad 
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Figure 7. (Color online) Transmission amplitude |T(Vg)| and phase a(V g ) for U/Y = 1 
and N = 4 equidistant levels with spacing S e . Decreasing 5 e /Y leads to a crossover 
from mesoscopic to universal behavior. 



level close to the chemical potential. The combination of these effects induces well 
separated transmission peaks, absent without interaction. The "universal" phase lapse 
can be understood as a result of the Fano effect of the effective renormalised levels, 
where it phase lapses appear in resonance phenomena (see lower panels in figure [7]). 
The mechanism leading to the phase lapses reflects also in the dot level occupancies, 
the broad level being filled and emptied via the narrow levels [14411145] . The detailed 
understanding of the charging of a narrow level coupled to a wide one in presence of 
interactions is of interest also in connection with experiments on charge sensing |146] . 
From an fRG analysis emerges a single-parameter scaling behavior of the characteristic 
energy scale for the charging of the narrow level with an interaction-dependent scaling 
function |147| . 



4.3. Outlook 

The described fRG approach presents a reliable and promising tool for the investigation 
of correlation and quantum fluctuation effects in quantum dots, the computational effort 
being comparable to a mean-field calculation. In contrast to the latter the fRG does 
not lead to unphysical artefacts. Besides the Kondo effect and transmission phases in 
multi-level dots, as explained above, also the competition between the Kondo effect 
and interference phenomena in more complex quantum-dot systems can be described 
J138pi48til51| . Furthermore, correlation effects for a large number of single-particle levels 
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can be discussed up to the limit of one-dimensional quantum wires. These applications 
involve inhomogeneous Luttinger liquids |152H155] , which will be discussed in section 
El as well as the influence of Luttinger liquid leads coupled to a quantum dot [156] . 
Challenging extensions and theoretical developments of the fRG involve the analysis of 
inelastic processes to describe dynamical properties, and the inclusion of higher-order 
contributions to access the physics at large U/Y. A fundamental question concerns the 
combination with the Keldysh formalism to describe nonequilibrium problems addressed 
in the next section. 

5. Quantum fluctuations in nonequilibrium and time evolution 

In this section we discuss quantum fluctuations in quantum dots or molecules in the 
presence of a finite bias voltage together with the time evolution into the stationary 
state. As in the previous section the quantum fluctuations are induced by the coupling 
to the reservoirs. If the maximum of temperature T and the distance to resonances 
S decreases, correlations effects from the Coulomb interaction lead to an increased 
renormalisation of the couplings and the excitation energies hi of the quantum dot 
(like e.g. magnetic field, level spacing, single-particle energy, etc.). Renormalisation 
group (RG) methods, which expand systematically in the reservoir-system coupling (in 
contrast to the fRG method described in the previous section HJ where an expansion in 
the Coulomb interaction is performed), show that these renormalisations are typically 
logarithmic or power laws. In nonequilibrium or for the time evolution, the voltage 
V and the inverse time scale 1/t are two new energy scales, which can cut off the 
RG flow. New physical phenomena emerge, which we will illustrate in this section by 
considering very basic two-level quantum dots, like two spin states (Kondo model) or 
two charge states (interacting resonant level model (IRLM)). Conceptually, one has to 
distinguish between the two different regimes of strong and weak quantum fluctuations. 
In strong coupling, an expansion in the renormalised coupling is no longer possible, 
and, except for the case of strong charge fluctuations in the IRLM or for the case 
of moderate Coulomb interactions (which can be treated with fRG methods), the 
nonequilibrium case remains a fundamental yet unsolved problem. Especially for the 
case of strong spin fluctuations, represented by the nonequilibrium Kondo model or the 
nonequilibrium single-impurity Anderson model, an analytical or numerical solution is 
still lacking. Even very basic questions, like e.g. the splitting of the Kondo peak in 
the spectral density by a finite bias voltage are not yet clarified. In weak coupling, 
where a controlled expansion in the renormalised coupling is possible, many results 
are already known. The simplest ansatz is to take the bare perturbation theory, as 
explained in section [21 use a Lorentzian broadening of the energy conservation (induced 
by relaxation and decoherence rates), and replace the bare couplings and excitation 
energies by the renormalised ones obtained from standard equilibrium poor man scaling 
(PMS) RG equations cut off by the maximum A c = max{T, hi, V, 1/t} of all physical 
energy scales. However, it turns out that this approach is not sufficient. In contrast 
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to temperature, the bias voltage is not an infrared cutoff since it tunes the distance to 
resonances, where new physical processes, like e.g. cotunneling (see section [2]), can occur. 
Therefore, the correct cutoff is not the voltage but the distance Si to resonances, where 
RG enhanced contributions occur. Furthermore, at resonance Si = 0, these contributions 
are cut off by relaxation and decoherence rates Tj, i.e. one should use the energy scale 
\Si + iTi\ as cutoff parameter for the couplings. An RG approach has to be developed to 
reveal how the various cutoff parameters influence the couplings. It turns out that the 
rates Tj are transport rates, which have to be determined from a kinetic (or quantum 
Boltzmann) equation, in contrast to rates describing the decay of a local wave function 
into a continuum. Therefore, the RG has to be combined with kinetic equations. The 
transport rates Tj depend also on voltage and are important parameters to prevent the 
system from approaching the strong coupling regime for voltages larger than the strong 
coupling scale T c . Furthermore, it turns out that terms can occur in the renormalised 
perturbation theory which are not present in the bare one, like e.g. the renormalisation 
of the magnetic field in linear order in the coupling for the Kondo model (similar effects 
can also happen in higher orders for generic models). For the time evolution, these 
terms are of particular interest in the short-time limit, since in this regime they are cut 
off by the inverse time scale 1/t and lead to universal time evolution. At large times, 
non-Markovian parts of the dissipative kernel of the kinetic equation lead to many other 
interesting effects. Among them are unexpected oscillation frequencies involving the 
voltage, unexpected decay rates, power-law behavior, and a different cutoff behavior at 
resonances compared to stationary quantities. In particular, for metallic reservoirs with 
a constant density of states, it can be shown that all local physical observables decay 
exponentially accompanied possibly by power-law behavior. Since all these effects occur 
already for very basic two-level systems, it has to be expected in the future that many 
more interesting effects will be found with respect to the physics of quantum fluctuations 
in nonequilibrium. 

From a technical point of view the description of correlated quantum dots in the 
presence of quantum fluctuations in nonequilibrium is a very challenging problem and a 
playground for the development of new analytical and numerical methods. Concerning 
numerical methods, numerical renormalisation group (NRG) with scattering waves 
[157] , time- dependent NRG (TD-NRG) [68lll58l4l60j . time- dependent density matrix 
renormalisation group (TD-DMRG) |161U161l4165j . quantum Monte Carlo (QMC) with 
complex chemical potentials |166j . QMC in nonequilibrium [167H169] . and iterative 
path-integral approaches (ISPI) |170j have been developed. However, the efficiency 
of these methods is often not satisfactory in the regime of either strong Coulomb 
interaction, large bias voltage, or long times. Exact analytic solutions exist for some 
special cases |162[ 117114174] and scattering Bethe-ansatz methods have been applied 
to the IRLM |175j . Other analytical methods are mainly based on RG approaches. 
Besides conventional PMS [176] . improved frequency-dependent RG schemes |177H179| . 
flow equations [180] . and fRG methods jl41U142l[T8TT4184j have been used. Whereas 
the latter expands systematically in the Coulomb interaction, we will introduce in this 
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section a formally exact RG method, which expands systematically in the reservoir- 
dot coupling. It uses the perturbation theory described in section [2] and allows for a 
direct calculation of the kernel of the kinetic equation fllj) in Liouville space, the so- 
called real-time RG method (RTRG) |185H187] , which recently has been technically 
optimized and formulated in pure frequency space (RTRG-FS) [5]. The method has 
been successfully applied to calculate analytically all static and dynamical aspects of 
weak spin fluctuations for the anisotropic Kondo model at finite magnetic field |188H190] 
and strong charge fluctuations for the IRLM [191] . Many aspects of these models can 
also be derived with other methods, which will be mentioned in the following at the 
appropriate places. 

5.1. Real-time RG in frequency space (RTRG-FS) 

The aim of RTRG-FS is to combine the diagrammatic expansion in Liouville space, 
as described in section |2} with RG to resum systematically infrared divergences. We 
consider a time translational invariant system and use as an initial condition at to = 
that the dot and the reservoirs are decoupled. The kinetic equation (jlj) can be written 
in Laplace space as 

* ) = 7^IfW p(t = 0) ' (16) 

where p(z) = J °° dt e lzt p(t) is the reduced density matrix of the dot in Laplace space 
and Lff (z) = Ld + S(z), where S(z) denotes the Laplace transform of the kernel 
E(t — t') = X(£, £'). The effective dot Liouville operator L^(z) contains all reservoir 
degrees of freedom and is of dissipative nature. The qualitative dynamics can be 
analyzed from the analytic structure of p(z), which is an analytic function in the upper 
half of the complex plane with poles and branch cuts only in the lower half. The 
single poles located at z l p = hi — iVi correspond to exponential decay with oscillation 
frequencies hi and decay rates Tj. From the pole at z — 0, the stationary state can 
be obtained via the solution of L^(i0 + )p = 0. Due to non-Markovian terms arising 
from the ^-dependence of the Liouvillian, the weight of these poles is changed by Z- 
factors and branch cuts can occur. For a constant density of states in the leads (normal 
metallic case) it can be shown generically |190U191j that the Z-factor leads to universal 
short-time behavior, whereas the branch cuts give rise to an exponential decay ~ e~ lz ^ 1 
accompanied by power-law behavior. z l h denotes the position of the branching point, 
which is shifted from some pole position by multiples of the electrochemical potentials of 
the reservoirs leading to oscillation frequencies involving the voltage. From the point of 
view of error correction schemes in quantum information processing it is quite important 
to understand these corrections to Markovian behavior |192H194] . 

It is unique to the RTRG-FS method that it provides direct access to the important 
quantity L£f (z). By systematically integrating out the energy scales of the reservoirs 
step by step, a formally exact RG equation can be derived for L^(z)\, where all reservoir 
energy scales beyond A are included. This RG equation is coupled to other RG equations 
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Figure 8. (color online) Two fundamental quantum dot models, (a) is the Kondo 
model, a spin-i coupled via exchange couplings J z ,± to two reservoirs, (b) is the IRLM, 
a spinless 1-level quantum dot coupled via tunneling rates T^r and Coulomb couplings 
U~l.r to two reservoirs. The electrochemical potentials are given by Ml/r = ±V/2 



for the couplings. Similar schemes can be developed for the calculation of the transport 
current [5] and correlation functions [189]. All RG equations involve resolvents similar 
to the one occurring in (Tl6l) where z is replaced by A together with other physical energy 
scales. As a consequence, it can be shown that, besides temperature, each term of the 
RG equation has a specific cutoff scale Aj , which is generically of the form 

Aj = \E + y^ rijUg. - hi + iTi\ = \5i +iVi\ . (17) 

3 

Here, E is the real part of the Laplace variable, rij are integer numbers, and \i a denotes 
the electrochemical potential of reservoir a. It shows that the cutoff scale is given by the 
distance Si to resonances. Furthermore, it provides the generic proof that, at resonance 
Si — 0, the cutoff scale is given by the corresponding rate IV This issue was under 
debate for some time because it was speculated that electrons tunneling in and out via 
the same reservoir correspond to low-energy processes, which could possibly lead to a 
strong coupling fixed point even in the presence of a finite bias voltage [195] . However, 
it was argued that volt age- induced decay rates prevent the system from approaching 
the strong coupling regime |177U196U197j . The microscopic inclusion of decay rates as 
cutoff scales into nonequilibrium RG methods was achieved within RTRG |185H187] . 
flow equation methods [180], and RTRG-FS [5]. 



5.2. Applications 

The two models used to illustrate the basic physics of spin and charge fluctuations are 
sketched in figure [HJ One model is the Kondo model at finite magnetic field h already 
discussed in section HI where a spin- 1/2 couples via anisotropic exchange couplings J z /± 
to the spins of two reservoirs. We have assumed a symmetric coupling to the leads and 
note that during the exchange it is also allowed that a particle is transferred between the 
reservoirs. The model results from the Coulomb blockade regime of a quantum dot with 
one level, where charge fluctuations are frozen out and only the spin can fluctuate. This 
leads to an effective band width of the reservoirs of the order of the charging energy U. 
Anisotropic exchange couplings can be realized for a molecular magnet, see section [2j 
The other model is the IRLM, where the quantum dot consists of a single spinless energy 
level at position e. The dot interacts with the reservoirs via tunneling processes, which 
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are characterized by the tunneling rates T a = 2-irp 0i \T ,\ 2 , with a = L, R. In addition, 
there is a Coulomb interaction u a between the first site of the reservoir leads and the 
quantum dot, which are characterized by the dimensionless parameter U a = p a u a . In 
the following we denote the bare parameters by a super-index (°). 

The two models have in common that due to spin/charge conservation the effective 
Liouvillian has the same matrix structure. There are three nonzero poles of the resolvent 
ffT6|) at Zp = —iY\ and z^ = ±h — iT 2 . For the Kondo model, Yi/ 2 corresponds to the 
spin relaxation/decoherence rate and h is the renormalised magnetic field. For the 
IRLM, Ti is the charge relaxation rate, T 2 describes the broadening of the local level, 
and h = e is the renormalised level position. Denoting the two eigenstates of the dot 
by ± =t, 1= 1,0, the matrix elements Lff (i0 + ) SStS > s > = — iss'W~ s of the Liouvillian 
involve the rates W s for the process — s — > s. The stationary occupations are given 
by p s = W s /W, which contain already most of the interesting nonequilibrium physics. 
Similar rates can be defined to calculate the current. We consider temperature T = 
from now on to reveal the physics of quantum fluctuations. 

5.2.1. Kondo model. We consider first the stationary case in the weak coupling regime 
A c = max{V, h} ^> Tk, where the Kondo temperature 7k = T c is the strong coupling 
scale for spin fluctuations. The renormalised couplings from PMS cut off at the 
scale A c are denoted by J z /±. For the isotropic case, they are explicitly given by 
J = l/(21n(A c /T K )), with T K ~ D^W^ ' e -V(u (0) ) (£) ~ U denotes the effective band 
width of the reservoirs) . 

In lowest order in J z /±, one obtains a "golden rule" like expression for the rates |188j 

W s = ^Jl{4h5.- + (V-h)6r a (V-h)} , (18) 

where h, V > and 6r(w) = 1/2 + (1/vr) arctan(o;/r) is a step function broadened by 
T. The latter corresponds to a Lorentzian broadening of the energy conservation law 
by quantum fluctuations. The spin relaxation/decoherence rates and the renormalised 
magnetic field are given by 

ri = |j£{4/i+(V-/O0r 12 (V-/O} , (19a) 

r 2 = ^r 1 + ^ 2 , (196) 

where Ti 2 = T\ — T 2 and 

h = (1 - J z + Jf) h<® 

-hJl\n-^—+ 1 -(V-h)Jlln A - , (20) 

\h + iT 12 \ 2 \V - h + iT 12 \ 

where unimportant terms ~ 0(J 2 ) without a logarithm have been left out for h. 

The expression for the rates W s up to 0(J 2 ) can be interpreted as follows. At 
high energies A > A c , the voltage and magnetic field are not relevant and the exchange 
couplings are renormalised according to the PMS equations. This leads to an effective 
band width A c for the reservoirs with renormalised exchange couplings cut off at A c . 
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Figure 9. g-factor g = 2dh/dhS°> for the isotropic Kondo model with V = 10 4 D, 



The precise value of A c is not important since the renormalisation is logarithmic. After 
this step one uses lowest order perturbation theory with broadened energy conservation 
leading to the result ( TT8l) for the rates. However, this interpretation does not work for 
the rates in 0(J 3 ) and it fails for the renormalised magnetic field h already in O(J) and 
0(J 2 ). We see that h contains a term linear in J z , which does not occur in perturbation 
theory, and was also discussed e.g. in |198| . It shows that it is generically not possible 
to calculate coefficients of certain orders in the renormalised couplings by comparing 
with bare perturbation theory. In addition, we find logarithmic corrections in 0(J 2 ). 
As was already discussed generically by Equation ( TTTj) . they occur at resonances h = 
or V — h. As was emphasized in |176j . they are of 0( J 2 In J) at resonance and remain a 
perturbative correction for J <C 1. They occur because the exchange couplings are not 
completely cut off by A c but partially by smaller cutoff scales. The calculation of the 
prefactor of the logarithmic terms is quite subtle and follows from the structure of the 
RG equations, where each term has its own cutoff scale Aj, see Equation ( JT7|) . Similar 
logarithmic corrections occur in 0(J 3 ) for the rates and for T1/2. In turns out that 
the logarithmic terms for h and Ti/ 2 contain an unexpected cutoff scale T^ = Ti — IV 
This is generic for all quantities entering the time evolution. It occurs because the pole 
positions of the resolvent ( TT6l) have to be calculated self-consistently. 

Logarithmic enhancements have been studied for the magnetic susceptibility and for 
the conductance also in |178[I179[I199] . using slave particle and flow equation methods. 
The logarithmic terms in the magnetic field have first been calculated in |188] using 
RTRG-FS. They lead to a suppression of the renormalised g-factor as function of h/V 
at h = and h — V, i.e. a nonequilibrium induced effect, see figure [9J Experimentally, 
it is proposed to be measured by using a three-lead setup, where the third lead is a 
weakly coupled probe lead |188| . A similar interesting nonequilibrium effect occurs 
for the logarithmic enhancement of the magnetic susceptibility x{h/V) at h = 0, as 
first proposed in |178j . All spin-spin correlation functions were calculated analytically 
in |189| using RTRG-FS. Here, the frequency variable of the correlation function enters 
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Equation (j!7p as well. At resonances, kink structures and logarithmic terms were found 
for the imaginary and real part of response functions, respectively. 

In |190j . the time evolution has been calculated from the inverse Laplace transform 
of Equation (TT6|) . In the short-time limit t <C A" 1 , the PMS equations are cut off by 
the energy scale 1/t leading to time dependent exchange couplings J^n- In O(J) only 
the Z-factor from the linear ^-dependence of L^(z) is important. In terms of the Pauli 
matrices cr, the local density matrix can be written as p(t) = \ + (S)(t) cr with 

(S)(t) = Z t (S)(0) = {l-2(Jt-j(V)}(S)(0) . (21) 

Inserting the solution of the PMS equations, one obtains universal logarithmic (power 
law) time evolution for the isotropic (anisotropic) case. A similar result has been found 
for the longitudinal spin dynamics in |200U201j for the special case of the ferromagnetic 
Kondo model, which was also confirmed by TD-NRG calculations [68] . In |190] it was 
also shown that the short-time behavior of the conductance can be calculated from the 
golden rule expression and replacing the exchange couplings by J*/ ± . 

In the long-time limit t ^> A" 1 , the cutoff scale for the PMS couplings is given by 
A c . In leading (Markovian) order one obtains the usual exponential behavior from the 
single poles of the resolvent f|T6|) . where the longitudinal/transverse spin decays with 
Ti/2 and the transverse spin oscillates with h. Interesting non-Markovian corrections 
occur in 0(J 2 ), where logarithmic contributions ~ J 2 {z — z l h ) ln(A c /(z — zfy) in Z^f (z) 
lead to branch cuts. The branching point of the logarithm is generically given by 
z l> = z i + n ^i l - e - * s shifted by multiples of the voltage from some pole position. As 
a consequence, an exponential behavior ~ J 2 e~ tZpt e~ mV is obtained, accompanied by 
power-law behavior ~ 1/t 2 from the branch cut integral. The result explains why the 
voltage occurs generically in the oscillation frequency, which is consistent with exact 
solutions at special Thoulouse points of two-level systems |173j . Furthermore, it shows 
that all terms decay exponentially but with unexpected decay rates and oscillation 
frequencies relative to the Markovian terms. This is illustrated in figure [10] by the 
time evolution of the transverse spin for the strongly anisotropic case T 2 ^> Ti, which 
is typical for molecular magnets. For short times, one obtains the expected result of 
oscillations with the magnetic field and decay with the spin decoherence rate T 2 . These 
terms decay quickly for large decoherence rate. For longer times, a crossover is obtained 
to an oscillation with the voltage V and a decay with the opposite spin relaxation rate Ti. 
Finally, we mention that close to resonances, i.e. for \S\ <C A c , 1/t, with \S\ = \V — h\,h, 
logarithmic terms ~ J 2 5tln \(5 + iTi2)t\ occur for the time evolution of the transverse 
spin. Similar to the logarithmic terms in h and Ti/ 2 , they are cut off by T 12 = Y\ — T 2 . 

5.2.2. Interacting resonant level model For the IRLM the PMS equations lead to a 
renormalised tunneling rate, given by the power- law T a = T a (D/A) 9a , with exponent 
g a = 2U a — ^o U 2 . Cutting off this scale at T c = T L + T R , one defines self-consistently 
the strong coupling scale T c for the importance of charge fluctuations. The level position 
e = e^ ) and the Coulomb interaction U a = Ua remain unrenormalised |202j . In lowest 
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Figure 10. (Color online) |Sa;(£)| = (|<Sa;(t)|) in the anisotropic Kondo model for 
V = 2h (0 *> = 100 T K and various values of the anisotropy c 2 = (J { z 0) ) 2 - (j| 0) ) 2 , with 
S x (0) = 1/2, 5*1/(0) = 0. The dips have their origin in the oscillations of S x (t). 



order in T a , one obtains a golden-rule like expression for the rates 



(22) 



where /xl/r = ±V^/2 are the electrochemical potentials of the leads. The energy 
broadening is given by 1^ = Ti/2 with the charge relaxation rate I\ = Yl a T«- 
Since a single tunneling process changes the charge, there is a unique cutoff parameter 
A° = |/i a — e + iT 2 \ for the PMS tunneling rates. This cutoff scale will be used in 
the following to define T a . For \/j, a — e| ~ T c , we are in the regime of strong charge 
fluctuations. We note the essential difference to the Kondo model, where several 
cutoff parameters are relevant. The current in lead a can be calculated from the rate 
equation I a = j-2n(W^po — W~p±) with the occupation probabilities p\ = W + /Y\ and 
Po = W~ /T\. These results are quite specific to the IRLM and are due to its elementary 
form. They have been conjectured in [2U2j and later were confirmed microscopically 
by RTRG-FS [191] and fRG [184] . It is surprising that these results are consistent 
with NRG [U2] and TD-DMRG [T62] calculations even for values of U a ~ 0(1), 
and capture all features obtained by field theoretical |162[ll74j and scattering Bethe 
ansatz [175] approaches. In particular, negative differential conductance is obtained for 
large voltages due to the power-law suppression of the tunneling rates. 

Defining A c = max{A^,A^} and expanding in U a , one finds logarithmic 
enhancements close to the resonances A" = 

9a ' A r 



'£ 



2U n In 



+ 0(U 2 



(23) 
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It is important to notice that if the level is in resonance with one of the reservoirs 
it is not with the other, i.e. exactly at resonance the cutoff scale is T 2 for one rate 
and the voltage V for the other rate. This has to be contrasted to speculations that 
both rates are cut off by the geometric average y/T 2 V |203j . which leads to incorrect 
power-law exponent for the on-resonance current as function of the voltage. Therefore, a 
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microscopic determination of the precise cutoff scales from nonequilibrium RG methods 
is very essential. 

The time evolution into the stationary state has been calculated for the IRLM 
in |191| . The basic features already found for the Kondo model were recovered, showing 
that these effects even hold for models with charge fluctuations. It turns out that the 
oscillation frequencies are generically given by hi + J^ ■ n^\i ap with integer numbers rij, 
i.e. the excitation energies h t of the dot are shifted by arbitrary multiples of the chemical 
potentials of the reservoirs. For the IRLM it turns out that the algebraic part of the 
time-evolution is ~ (l/t) 1_9a , i.e. the exponent depends on the Coulomb interaction 
strength. 

5.3. Outlook 

The status of the field of quantum fluctuations in nonequilibrium for strongly correlated 
quantum dots is that powerful techniques have been developed to study the weak- 
coupling limit in a controlled way. Elementary models of spin and charge fluctuations 
are well understood in this limit and the methods can now be applied to more complex 
quantum dot models, such as discussed in section 12.21 However, two fundamental 
issues are still open. First, the case of strong spin fluctuations, represented in its most 
elementary form by the isotropic Kondo model in the limit max{T, h, V} ~ Tk, is one 
of the most fundamental unsolved problems. The case of strong charge fluctuations 
at resonances has so far only been understood for the IRLM, where it seems that 
the broadening of the level together with a PMS renormalisation of the tunneling 
couplings captures the essential physics. Whether this holds also for more complicated 
models has to be tested in the future. Secondly, most of the methods used to 
describe nonequilibrium properties of quantum dots are parametrized by the many-body 
eigenstates of the isolated quantum dot. Therefore, they can not be extended easily to 
larger systems like multi-level quantum dots or quantum wires. Two exceptions are 
the flow-equation method [180j and nonequilibrium fRG methods. |141U142[ll81H184j . 
In fRG, as already described in section HI a perturbative expansion in the Coulomb 
interaction is used and the vertices are parametrized in terms of single-particle levels. 
Therefore, although fRG is restricted to the regime of moderate Coulomb interactions, 
the potential of the method lies particularly in the possibility to treat multi-level 
quantum dots and quantum wires. Preliminary studies have tested the fRG for the IRLM 
and the single-impurity Anderson model (SIAM). For the IRLM, a static version was 
used, where the frequency dependence of the Coulomb vertex was neglected |184U191] . 
An agreement was found with RTRG-FS and TD-DMRG |162j . Recently, a dynamic 
scheme with frequency-dependent 2-particle vertices has been developed and was used 
to analyze the nonequilibrium SIAM in the regime of strong spin fluctuations [14111142] . 
A good agreement with TD-DMRG, ISPI and QMC results was obtained for moderate 
Coulomb interactions |204j . As will be described in the next section, the static version 
of nonequilibrium fRG has also been applied to quantum wires and it is a challenge for 
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future research to generalize the dynamic fRG scheme to quantum wires as well. 

6. Correlation effects in quantum wires 

6.1. An introduction to Luttinger liquid physics 

A different type of correlation physics than discussed so far is found in (quasi) one- 
dimensional (Id) quantum wires at low temperatures. By this we mean electron systems 
confined in two spatial directions such that only the lowest Id subband is occupied. 
Compared to the quantum dots considered above they contain many correlated degrees 
of freedom and the single-particle level spacing of the wire becomes the smallest energy 
scale ("almost" continuous spectrum; see below). The dominant quantum fluctuations 
are now driven by the interaction U itself. In fact, below we will first discuss isolated, 
translationally invariant wires. Still, as discussed in the main part of this section, 
the physics becomes even more interesting if the coupling to leads (or the coupling 
between two wires) is included. Quantum wires are realized in single-walled carbon 
nanotubes [2051 I206J . specifically designed semiconductor heterostructures |207U208j . 
and in atom chains which form on certain surfaces |209j . A different class of quasi 
Id systems are highly anisotropic (chain-like) bulk materials (for a recent review see 
e.g. |210j ). Compared to the first type of systems it is less clear how a single wire of 
this class can be incorporated as an element in an electronic transport nanodevice and 
they will thus play a minor role in the present discussion. 

Although being far from experimentally realizable, translationally invariant models 
of Id electrons with sizeable two-particle interaction (Coulomb repulsion) were already 
studied in the fifties and sixties. Exactly solvable models were constructed leading 
to a rather quick (compared to higher- dimensional correlated systems) gain in the 
understanding of correlation effects. Tomonaga showed that all excitations of a Id 
correlated electron system at (asymptotically) low energies are collective in nature 
("plasmons") rather than single-particle-like as in the three dimensional counterpart 
(Landau quasi-particles) |211j . In particular, collective spin and charge density 
excitations travel with different speed. Quite often this difference in velocities of 
collective excitations is mistaken as being characteristic for Id systems although it might 
also occur in higher dimensional Fermi liquids. Only when adding, that no additional 
quasi-particle excitations are possible spin-charge separation becomes a unique feature 
of Id correlated electrons. As a second (related) interaction effect certain correlation 
functions display power-law behavior as a function of energy u and momentum k as first 
shown for the single-particle momentum distribution function n(k) by Luttinger [212] 
and Mattis and Lieb |213j . Using the concepts of statistical physics one can view Id 
chains of electrons as being (quantum) critical, although the exponents depend on the 
details of the underlying single-particle model (band structure, band filling) as well as on 
the interaction strength. Solyom [214] and Haldane [215j showed using renormalisation 
group (RG) arguments that the physics of the two exactly solvable models — being 
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quite similar they are nowadays considered as a single model, the Tomonaga-Luttinger 
model — is generic for Id interacting electrons as long as correlations do not drive the 
system out of the metallic phase (e.g. into a Mott insulator phase). This led to the term 
Luttinger liquid (LL) physics for the above correlation effects. In particular, Haldane 
argued that all the exponents appearing in correlation functions of a spin-rotational 
invariant system with repulsive two-particle interaction can be expressed in terms of 
a single model parameter dependent number K < 1, with K — 1 for noninteracting 
electrons. Together with the velocities of the charge- and spin-density excitations v c and 
v 8 , K completely determines the low-energy physics. For a given microscopic model the 
strategy to obtain the low-energy physics is thus the following: one has to determine 
the three Luttinger liquid parameters as functions of the model parameters (ways to 
achieve this are e.g. discussed in |216j ) and plug them into the analytic expressions for 
correlation functions computed within the Tomonaga-Luttinger model. 

Despite the intense effort (see e.g. [205H210] ) a commonly accepted experiment 
which shows LL behavior beyond any doubts is still pending. This has to be contrasted 
to the situation in quantum dots discussed earlier where the appearance of Kondo physics 
as a consequence of local correlations has convincingly been shown (see above). The 
lack of clear cut experiments is partly related to the status of the theory. LL theory 
only makes predictions for the asymptotic low-energy behavior but does not provide 
the scales on which this sets in. In the construction of the Tomonaga-Luttinger model 
all microscopic scales as e.g. given by the band curvature |217] and the shape of the 
two-particle interaction |218j are disregarded. One can imagine that the upper energy 
scale T c beyond which LL physics can be found becomes so small that power laws 
are masked by other scales such as the single-particle level spacing set by the length 
of the wire. A profound comparison to experiments requires a knowledge of T c and 
possible lower bounds for LL physics while methods which allow to extract these scale 
for microscopic models are rare. It is thus mandatory to develop new methods which 
capture LL physics and can directly be applied to microscopic models. Further down, 
we will return to this issue. We note in passing that the standard ab initio method — 
density functional theory — cannot be applied as the existing approximation schemes fail 
when it comes to LL physics. 

When using a quantum wire as part of a low-temperature nanodevice one has to 
be alert that correlation effects might alter the electronic properties. To understand 
this let us first imagine an electron tunneling from the tip of a scanning tunneling 
microscope into the bulk of a Id quantum wire. Using Golden Rule-like arguments it 
becomes clear that the differential conductance dl/dV computed from the tunneling 
current / is determined by the product (more precisely the convolution) of the single- 
particle spectral function p(cu) (the "local density of states") of the wire and the tip. 
As the latter is constant at small energies, dl/dV displays the LL power law of the 
wire spectral function p(u) ~ \u — fi\ ahulk , with the wires chemical potential /i and 
«buik = {K + K" 1 — 2)/4 |216| . The exponent changes to a en d = (K _1 — l)/2 if the 
electron tunnels into the end of a LL as it is e.g. the case in end-contacted wires |216j . 
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Evidently, tunneling in and out of either the bulk or the end of a wire plays a crucial 
role when using it in a nano device. 

Secondly, we imagine a local inhomogeneity in the quantum wire. It might either 
appear (in a difficult to control way) in the process of producing the wire or might 
intentionally (in a controlled way) be introduced to design a functional device. As 
examples for the latter one can think of double barriers |219] defining a dot region with 
energy levels that can be tuned via a backgate (see the discussion on transport through 
quantum dots) or junctions of several quantum wires |220U221j . One can easily imagine 
that because of the collective nature of the excitations in a LL any inhomogeneity with 
a (single-particle) 2/cF-backscattering component, with fcp being the Fermi momentum, 
strongly affects the low-energy physics. A first step to substantiate this expectation is to 
compute the 2fcp-component of the static density-density-response function x(o) of the 
Tomonaga-Luttinger model [222] . For a noninteracting Id electron system it diverges 
logarithmically when q approaches 2kp (Lindhard function). The divergence is enhanced 
and becomes power-law like x(l) ~ \Q ~ 2/cf| a ' _1 if the interaction is turned on. This 
shows that the inhomogeneity strongly couples to the system and linear response theory 
breaks down. This physics was further investigated applying different RG methods to 
different models [223H227] . These studies show that even a single impurity at low 
energy scales effectively acts as if the chain is cut at the position of the impurity 
with open boundary conditions at the two end points. E.g. for temperatures T — > 
the linear conductance G(T) across the impurity is suppressed in a power-law fashion 
G(T) ~ T 2oicnd as can again be understood using Fermis Golden Rule and the power- 
law scaling of the local spectral function at the end of a quantum wire (tunneling 
from end to end of two LLs). Remarkably, the exponent is independent of the bare 
scattering potential. Also a Breit-Wigner resonance of G as a function of the level 
position (modified by) V^ showing up in the double-barrier geometry is strongly altered 
by the interaction. In case of a "perfect" resonance with peak conductance e 2 /h (per 
spin), that is for equal left and right barriers (which is experimentally difficult to realize), 
the resonance width becomes zero, while the resonance completely disappears, that is 
G — > for all V s , in all other cases. 

The discussed physics of inhomogeneous LLs can be understood in terms of an 
effective single particle problem. Using RG techniques one can show that the interplay 
of a local inhomogeneity and two-particle correlations leads to an effective oscillatory 
and slowly decaying (~ 1/x, with the distance x from the inhomogeneity) single-particle 
scattering potential of range 1/6, where 8 is the largest of the relevant energy scales 
(e.g. temperature). The wave-length of the potential is set by the (common) chemical 
potential /x of the leads. Scattering off this (so-called Wigner-von Neumann potential) 
leads to the discussed power laws in transport with the exponent given by the amplitude 
of the potential. [221H22T] 

In the remaining part of this section on coherent transport through quantum wires 
we describe four examples of recent attempts (by three of the present authors together 
with varying colleagues) to gain a detailed understanding of the interplay of correlations 
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Figure 11. (Color online) Temperature dependence of the linear conductance for the 
extended Hubbard model with 10 4 sites and a single site impurity ( "6- impurity" ) of 
strength V = 10, for a Hubbard interaction U = 1 and various choices of U'\ the 
density is n = 1/2, except for the lowest curve, which has been obtained for n = 3/4 
and U' = 0.65 (leading to a very small backscattering interaction); the dashed line is 
a power law fit for the latter parameter set. 



and local inhomogeneities in systems with different geometries. We believe that our 
insights should be kept in mind when designing transport setups to search for LL physics. 
They might also be of relevance when quantum wires (in the above defined sense) are 
used as elements in future nanoelectronic devices. In this case correlation effects might 
either be used to enhance the functionality or must be tuned away if they corrupt the 
latter. 

In all the examples we use the functional renormalisation group (fRG) method [137] 
(for a brief review see also [228] ) to (approximately) treat the Coulomb interaction. It 
was already introduced in section [U Compared to other methods it has the advantages 
that it (i) can directly be applied to microscopic models (continuum or lattice; flexible 
modeling), (ii) captures all energy scales (not restricted to the asymptotic low-energy 
regime), and (iii) can be applied to systems with many correlated degrees of freedom. 
In the present context (iii) means chains of realistic length (in the micrometer range) 
coupled to effectively noninteracting (Fermi liquid) leads. The aspects (i) and (ii) 
together assure that T c and additional energy scales affecting LL physics can be 
investigated. As mentioned in section [H the approximations inherent to the fRG 
approach are justified up to moderate renormalised two-particle interactions. In the 
implementation (static self-energy) used for quantum wires (i) the LL exponents come 
out correctly (only) to leading order in the two-particle interaction and (ii) inelastic 
processes generated by the two-particle interaction are neglected. We will return to the 
latter in the outlook. 
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6.2. The role of two-particle backs cattering 

In Id systems at low energies and incommensurate fillings only forward (momentum 
transfer q rs 0) and backward (momentum transfer q ps 2kp) two-particle scattering is 
active. All other processes are suppressed due to momentum conservation and phase 
space restrictions (Fermi surface consists of only two points). The relative importance of 
the two processes is determined by the (real-space) range of the two-particle interaction, 
which in turn is set by the screening properties of the "environment" , e.g. the substrate 
on which the carbon nanotube is placed. If screening is strong, the two-particle 
interaction becomes short ranged in real space and backscattering is sizable and vice 
versa. Backscattering processes involving a spin up and a spin down electron are not 
part of the Tomonaga-Luttinger model. For a translationally invariant LL it was shown 
using an RG analysis |214] that these processes do not modify the asymptotic low- 
energy physics (they are "RG irrelevant"), but that they become irrelevant only on a 
scale which is exponentially small in the bare backscattering strength (they flow to zero 
only logarithmically). From this one can expect that the LL power laws only occur on 
exponentially small scales T c also for inhomogeneous wires. For the Hubbard model with 
open boundaries it was shown that T c ~ exp (—ttvf/U), where U is the local Coulomb 
repulsion and v-p the Fermi velocity |229j . In figure [11] this behavior is exemplified for 
the T-dependence of the linear conductance G of a wire with a single impurity described 
by the extended Hubbard model [155j . This lattice model consists of a standard tight- 
binding chain with nearest-neighbor hopping t, a local two-particle interaction U as 
well as a nearest-neighbor one U' . Throughout the rest of this section we use t as the 
unit of energy (that is we set t = 1). For a fixed band filling the relative strength 
of the forward and backward scattering can be modified by varying the ratio U'/U. 
The considered size of iV = 10 4 lattice sites corresponds to wires in the micrometer 
range, which is the typical size of quantum wires available for transport experiments. 
For U' = (strong backscattering) due to logarithmic corrections the conductance 
increases as a function of decreasing T down to the lowest temperatures in the plot, 
T c being smaller than the latter. For increasing nearest-neighbor interactions U' the 
(relative) importance of backscattering decreases and a suppression of G(T) at low T 
becomes visible, but in all the data obtained at quarter-filling n — 1/2 the suppression 
is much less pronounced than what one expects from the asymptotic power law with 
exponent 2a en d. By contrast, the suppression is much stronger and follows the expected 
power law more closely if parameters are chosen such that two-particle backscattering 
becomes negligible at low T, as can be seen from the conductance curve for n = 3/4 and 
U' = 0.65 in figure [TTJ The value of K for these parameters almost coincides with the 
one for another parameter set in the plot, n = 1/2 and U' = 0.75, but the behavior of 
G(T) is completely different. To avoid interference effects discussed next, in the present 
setup the two leads are coupled "adiabatically" to the wire, such that in the absence 
of the impurity the conductance would be 2e 2 /h for all temperatures smaller than the 
bandwidth [227]. 
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One can conclude that for wires in the micrometer range clear power laws can 
only be observed if screening is weak, leading to small two-particle backscattering 
[1551I2241 - I226J . Note that at T ~ irvp/N finite size effects set in (the level structure 
becomes apparent), as can be seen at the low T end of some of the curves in the figure. 

6. 3. Overhanging parts 

In most transport experiments the wires are not end-contacted and the leads do not 
terminate at the contacts. This has to be contrasted to the modeling in which almost 
exclusively end-to-end contacted wires are considered. It is thus crucial to understand 
how overhanging parts of the wire and the leads alter the transport characteristics. 
Before tackling this problem we have to discuss transport through an interacting wire 
with two contacts — above we only mentioned the single contact case with power- 
law scaling G(T) ~ T acnd . One might be tempted to argue, that two contacts can 
be understood as two resistors in series, which then have to be added leading to 
G(T) ~ T acnd also for a (clean) wire with two contacts to semi-infinite leads. Although 
the result is correct, as can be shown using (coherent) scattering theory [22711230] . this 
simple argument ignores, that adding of the resistance (of resistors in series) only holds 
if the transport is incoherent, which is not the case at low temperatures as considered 
here. 

We here exemplify the role of overhanging parts by considering overhanging leads. 
In figure [12] G(T) is shown for a spinless tight-binding wire with nearest-neighbor 
interaction U' coupled to two semi-infinite leads via tunneling hoppings Tl/r. For 
spinless models the low-energy LL physics can be described by K and the charge 
velocity v c . Several correlation functions are again characterized by power-law scaling. 
Compared to models with spin the analytic dependence of exponents on K is modified 
(see e.g. [216] ). Three setups with different overhanging parts are considered (see the left 
inset). As becomes clear each non-zero number of overhanging lattice sites N lead infers a 
new energy scale kv-f/NJ^ and the simple power law (with exponent a e nd; dashed line) 
is divided in subsections separated by extensive crossover regimes. If clearly developed 
the power laws in the subsections all have exponent a en d (see the inset of figure [T2l . 
Overhanging parts of the wire lead to the same effect [231j . There is no reason to believe 
that adding the spin degree of freedom will change this. The same holds for the two 
effects discussed next and we thus stick to the spinless lattice model for the rest of this 
section. 

One can conclude that in a generic experimental setup with overhanging parts the 
power law is much more difficult to observe than suggested by most theoretical studies 
considering end-to-end coupling. 

6.4- A wire with two contacts and a single impurity 

We now exemplify the very rich interference effects which occur when considering 
coherent linear transport through interacting wires with several inhomogeneities. A 



Charge transport through single molecules, quantum dots, and quantum wires 



39 



0.2 



0.1 



o 



0.03 



110 



8- -0.2 - 



\ I 



V 



,20 1086 



Nio° 
\ 



\ ' ioVi 



10" 



10" 



10 
T 



10" 



10 



Figure 12. (Color online) Main plot: Linear conductance G of an interacting wire 
with U' = 0.5 of length TV = 10 4 + 1 as function of the temperature T. The coupling to 
the leads is located at the end of the wire and has a small amplitude Tr = XL = 0.25. 
Dashed line: Overhanging parts j% o ^ d = 0, that is end-to-end contacted wire. Dashed- 
dotted line: iV£ ad = 0,iV£ ad = 110. Dotted line: iV£ ad = 20,A/£ ad = 1086. The 
vertical dashed lines terminating at the different curves indicate the corresponding 
crossover scales. At the lowest scale T ~ 10~ 4 the power law is cut off by the level 
spacing ~ 1/iV of the disconnected wire. Left inset: Setups studied (excluding the 
end-to-end contacted one). Right inset: Effective exponents obtained from taking the 
logarithmic derivative of G(T). The solid horizontal line indicates a e nd- 



simple setup of this class is a wire with two end-to-end tunneling contacts to semi- 
infinite leads and an additional impurity in its bulk. As mentioned in the last subsection 
the concept of "adding resistances" cannot be applied in the present context. Using 
averaging over the (two) scattering phases one can show that for sufficiently large T (of 
the order 10 -2 in figure IT3]) the total linear conductance is given by G 



VG1G2G3 

with the conductances G?j of the individual barriers. For the present setup we have 
Gi/ 3 (T) ~ T a - d and G 2 (T) ~ T 2Q - d , leading to G(T) ~ T 2a - d |j. As is shown in 
figure [13] the dashed-dotted curve obtained by phase averaging nicely follows the bare 
data (solid and dashed line) down to T p w O(10~ 2 ) and a power law with the expected 
exponent develops for T p ;$ T ^C D. For T < T p details of the relative energy- level 
structure of the two "dots" defined by the three barriers matter and might even lead to 
resonances (see the solid line). The "dots" energy-level structure is set by the position 
of the bulk impurity and the solid and dashed line in figure [13] display data obtained 
for different positions. The example clearly shows that interference due to multiple 
inhomogeneities can set a lower bound to power-law scaling of the conductance. 



6. 5. Transport at finite bias voltages 

As our final example we study transport at finite bias voltages. In particular we 
consider the steady state current through a clean wire end-to-end coupled to two 



If Accidentally the exponent 2a onc i would also follow from the unjustified adding of resistances. 
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Figure 13. (Color online) Linear conductance G of an interacting tight-binding wire 
with tunnel couplings Tl = 2r = 0.1, an additional hopping impurity (single bond 
with reduced hopping) of strength if = 0.1, and nearest-neighbor interaction U' = 1 of 
length N = 10 4 as function of the temperature T. The solid and dashed curve show 
data obtained for different positions of the bulk impurity. Dashed-dotted line: result 
obtained by phase averaging. 



leads (reservoirs) which are hold on different electrochemical potentials /il/r = ±V/2. 
Applying a generalization of the fRG to non-equilibrium transport |183] one can 
show that in the limit of high tunnel barriers the effective single-particle potential 
generated by the interplay of the contacts and the interaction is a superposition 
of two decaying oscillations with wave-lengths set by /il and /x R . The respective 
amplitudes are proportional to Tl, R /(Tl + T R ) = IY/r^IY + Tr), with T a ~ Tj 
defined as in the preceding sections. Accordingly, the non-equilibrium local spectral 
function close to either of the two contacts displays power-law suppressions at /xl and 



Vte., pH 



\U) 



/iL/R.| aL/R , with exponents proportional to the respective couplings 



«l/r = rL/Rtt cn d/(rL + T R ). In contrast to the linear response regime the exponents 
now dependent on the strength of the inhomogeneities. This behavior is shown in 
figure [TH The spectral weight remains finite even at /il/r as the power-laws are cut 
off by the single-particle level spacing ~ 1/N of the disconnected wire. A transport 
geometry in which the two power laws show up was proposed in |183| . 



6. 6. Outlook 

This short review on correlation effects in Id quantum wires gives a flavor of the rich 
and partly surprising physics resulting from the interplay of the two-particle interaction 
and local (single-particle) inhomogeneities. An important ingredient missing in the 
current description using the (approximate) fRG for microscopic models are inelastic 
scattering processes resulting from the (Coulomb) interaction and leading to decoherence 
and dephasing at intermediate to large energies. We note that they are also only partly 
kept in the standard modeling which is based on the Tomonaga-Luttinger model. One 
can expect that inelastic processes set an upper energy scale T c for power-law scaling 
smaller than the maximal scale given by the band width. They will also play a prominent 
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Figure 14. (Color online) Solid line: Local spectral function (close to the left contact) 
in the steady-state of an interacting tight-binding wire of length TV — 2 • 10 4 with tunnel 
couplings Tl = 0.075, Tr = 0.15, and nearest-neighbor interaction U' = 0.5 as function 
of ui. The voltage is V = 0.1 



role when further considering non-equilibrium transport. First attempts to include 
such inelastic two-particle scattering for quantum dots in and out of equilibrium were 
mentioned in section H] and are discussed in [139H141] . 



List of symbols 

In this section we provide a list of the most relevant symbols used: 



Ae 

r 
v g 

V 

dl/dV 

t 

T K 

U 

h = k B 

T 

\s) 

H D 

Ht 

rpss' 

akcr 



a 
a - 



L,R 



"dako 
k 



level spacing 

tunnel coupling 

gate voltage 

bias voltage 

differential conductance 

time 

Kondo temperature 

charging energy 

units 

temperature 

exact many-body level spectrum and states 

dot Hamiltonian 

tunneling Hamiltonian 

tunnel amplitudes 

spin 

reservoirs 

reservoir Hamiltonian 

electron number 

label of the orbital state 
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Pa 


density of states 


Pa 


electro-chemical potential 


P (t) 


dot density operator 


L D 


dot Liouville operator 


£(*,*') 


transport kernel 


la 


tunnel current 


e* 


level energy 


N 


charge number 


u 


frequency 


D N 


magnetic parameter 


X 


electron vibration coupling 


n 


dwell time of electrons in the system 


Q 


pumped charge 


G(o) 


(bare) Green function of the dot 


T,(ioj) 


self-energy 


G 


conductance 


Si 


distance to resonances 


Zi 


poles of the reduced dot density matrix p in Laplace space 


h 


magnetic field 


Jz/± 


exchange couplings 


W s 


rates for the process —s^-s 


X 


susceptibility 


A 


flow parameter 


Ac 


cutoff 


K 


Luttinger liquid parameter 


v c / s 


velocities of the charge/spin-density excitations 


kp, vf 


Fermi momentum, Fermi velocity 


U' 


nearest-neighbor interaction 


ttbulk 


power-law exponent of the wire spectral function 


"end 


power-law exponent of an electron tunneling into the end of a LL 
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